General info

I created this website to track the progress with the sweetgum genetic study.

There are 7 parts:

  1. Libraries, dataset loading, and data wrangling.
  2. SAS code.
  3. Survival rate plots.
  4. Groundline diameter (GLD) plots and stats.
  5. Tree height plots and stats.
  6. DBH plots and stats.
  7. Live crown height plots and stats.

Data Wrangling

library(tidyverse)
## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr     1.1.4     ✔ readr     2.1.4
## ✔ forcats   1.0.0     ✔ stringr   1.5.1
## ✔ ggplot2   3.4.4     ✔ tibble    3.2.1
## ✔ lubridate 1.9.3     ✔ tidyr     1.3.0
## ✔ purrr     1.0.2     
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors
library(ggpubr)
library(plotly)
## 
## Attaching package: 'plotly'
## 
## The following object is masked from 'package:ggplot2':
## 
##     last_plot
## 
## The following object is masked from 'package:stats':
## 
##     filter
## 
## The following object is masked from 'package:graphics':
## 
##     layout
library(htmltools)

hane_crown <- read.csv("https://raw.githubusercontent.com/azd169/R-datasets/main/hane_crown.csv",
                       header = T)
hane_gld <- read.csv("https://raw.githubusercontent.com/azd169/R-datasets/main/hane_gld.csv",
                       header = T)
hane_dbh <- read.csv("https://raw.githubusercontent.com/azd169/R-datasets/main/hane_dbh.csv",
                       header = T)
hane_ht <- read.csv("https://raw.githubusercontent.com/azd169/R-datasets/main/hane_ht.csv",
                       header = T)
hane_survival_year <- read.csv("https://raw.githubusercontent.com/azd169/R-datasets/main/hane_survival_year.csv",
                       header = T)
hane_survival_year_family <- read.csv("https://raw.githubusercontent.com/azd169/R-datasets/main/hane_survival_year_family.csv",
                       header = T)
hane_survival_year_location <- read.csv("https://raw.githubusercontent.com/azd169/R-datasets/main/hane_survival_year_location.csv",
                       header = T)
hane_survival_year_location_family <- read.csv("https://raw.githubusercontent.com/azd169/R-datasets/main/hane_survival_year_location_family.csv",
                       header = T)

hane_crown <- hane_crown %>%
  mutate(Year = fct_relevel(Year, c("Fifth", "Sixth", "Seventh")))

hane_gld <- hane_gld %>%
  mutate(Year = fct_relevel(Year, c("Initial", "First", "Second")))

hane_dbh <- hane_dbh %>%
  mutate(Year = fct_relevel(Year, c("Third", "Fourth", "Fifth", "Sixth", "Seventh")))

hane_ht <- hane_ht %>%
  mutate(Year = fct_relevel(Year, c("Initial", "First", "Second", "Third", "Fourth", "Fifth", "Sixth", "Seventh")))

hane_survival_year <- hane_survival_year %>%
  mutate(Year = fct_relevel(Year, c("Initial", "First", "Second")))

hane_survival_year_family <- hane_survival_year_family %>%
  mutate(Year = fct_relevel(Year, c("Initial", "First", "Second")))

hane_survival_year_location <- hane_survival_year_location %>%
  mutate(Year = fct_relevel(Year, c("Initial", "First", "Second")))

hane_survival_year_location_family <- hane_survival_year_location_family %>%
  dplyr::mutate(ID = paste(Location, Year, sep = ", ")) %>%
  mutate(ID = fct_relevel(ID, c("Hill Farm, Initial", "Hill Farm, First", "Hill Farm, Second",
                                "LA Tech, Initial", "LA Tech, First", "LA Tech, Second")))

SAS Code

Procedure:

External trees were removed from the analysis and all the Variables were averaged at plot level. A unique ID was created combining Genotype, Location, and Block (GLD in the SAS code)

Datasets used:

  1. WORK.DBH for DBH
  2. WORK.GLD for GLD
  3. WORK.HEIGHT for tree height
  4. WORK.CROWN for crown height

Here are the SAS models:

  1. GLD

PROC GLIMMIX DATA = GLD H PLOTS = ALL; CLASS Year Location Block Family GLB; MODEL GLD = Year|Location|Family / SOLUTION DDFM = KR; RANDOM intercept / SOLUTION SUBJECT = GLB TYPE = AR(1) CL; LSMEANS Location Family Year * Location Year * Family/ ADJUST = TUKEY LINES; ODS GRAPHICS OFF; RUN;

  1. Height

PROC GLIMMIX DATA = HEIGHT H PLOTS = ALL; CLASS Year Location Block Family GLB; MODEL Height = Year|Location|Family / SOLUTION DDFM = KR; RANDOM intercept / SOLUTION SUBJECT = GLB TYPE = AR(1) CL; LSMEANS Family/ ADJUST = TUKEY LINES; SLICE Year * Location / SLICEBY = Year ADJUST = TUKEY LINES DIFF; SLICE Year * Family / SLICEBY = Year ADJUST = TUKEY LINES DIFF; ODS GRAPHICS OFF; RUN;

  1. DBH

PROC GLIMMIX DATA = DBH H PLOTS = ALL; CLASS Year Location Block Family GLB; MODEL DBH = Year|Location|Family / SOLUTION DDFM = KR; RANDOM intercept / SOLUTION SUBJECT = GLB TYPE = AR(1) CL; LSMEANS Family/ ADJUST = TUKEY LINES; SLICE Year * Location / SLICEBY = Year ADJUST = TUKEY LINES DIFF; SLICE Year * Family / SLICEBY = Year ADJUST = TUKEY LINES DIFF; ODS GRAPHICS OFF; RUN;

  1. Crown Height

PROC GLIMMIX DATA = CROWN H PLOTS = ALL; CLASS Year Location Block Family GLB; MODEL CROWN = Year|Location|Family / SOLUTION DDFM = KR; RANDOM intercept / SOLUTION SUBJECT = GLB TYPE = AR(1) CL; LSMEANS Family/ ADJUST = TUKEY LINES; SLICE Year * Location / SLICEBY = Year ADJUST = TUKEY LINES DIFF; SLICE Year * Family / SLICEBY = Year ADJUST = TUKEY LINES DIFF; ODS GRAPHICS OFF; RUN;

SAS Output GLD

htmltools::includeHTML("https://raw.githubusercontent.com/azd169/SAS_outputs/main/hane-results-GLD.html")
## Warning: `includeHTML()` was provided a `path` that appears to be a complete HTML document.
## ✖ Path: https://raw.githubusercontent.com/azd169/SAS_outputs/main/hane-results-GLD.html
## ℹ Use `tags$iframe()` to include an HTML document. You can either ensure `path` is accessible in your app or document (see e.g. `shiny::addResourcePath()`) and pass the relative path to the `src` argument. Or you can read the contents of `path` and pass the contents to `srcdoc`.
Results: hane.sas

Results: hane.sas

The GLIMMIX Procedure

The Glimmix Procedure

Model Information

Model Information
Data Set WORK.GLD
Response Variable GLD
Response Distribution Gaussian
Link Function Identity
Variance Function Default
Variance Matrix Blocked By GLB
Estimation Technique Restricted Maximum Likelihood
Degrees of Freedom Method Kenward-Roger
Fixed Effects SE Adjustment Kenward-Roger

Class Level Information

Class Level Information
Class Levels Values
Year 3 First Initial Second
Location 2 Hill Farm LA Tech
Block 4 1 2 3 4
Family 6 AGH2 AGH25 AGHS1 AGHS2 AGHS3 AGHS4
GLB 48 AGH25_Hill Farm_1 AGH25_Hill Farm_2 AGH25_Hill Farm_3 AGH25_Hill Farm_4 AGH25_LA Tech_1 AGH25_LA Tech_2 AGH25_LA Tech_3 AGH25_LA Tech_4 AGH2_Hill Farm_1 AGH2_Hill Farm_2 AGH2_Hill Farm_3 AGH2_Hill Farm_4 AGH2_LA Tech_1 AGH2_LA Tech_2 AGH2_LA Tech_3 AGH2_LA Tech_4 AGHS1_Hill Farm_1 AGHS1_Hill Farm_2 AGHS1_Hill Farm_3 AGHS1_Hill Farm_4 AGHS1_LA Tech_1 AGHS1_LA Tech_2 AGHS1_LA Tech_3 AGHS1_LA Tech_4 AGHS2_Hill Farm_1 AGHS2_Hill Farm_2 AGHS2_Hill Farm_3 AGHS2_Hill Farm_4 AGHS2_LA Tech_1 AGHS2_LA Tech_2 AGHS2_LA Tech_3 AGHS2_LA Tech_4 AGHS3_Hill Farm_1 AGHS3_Hill Farm_2 AGHS3_Hill Farm_3 AGHS3_Hill Farm_4 AGHS3_LA Tech_1 AGHS3_LA Tech_2 AGHS3_LA Tech_3 AGHS3_LA Tech_4 AGHS4_Hill Farm_1 AGHS4_Hill Farm_2 AGHS4_Hill Farm_3 AGHS4_Hill Farm_4 AGHS4_LA Tech_1 AGHS4_LA Tech_2 AGHS4_LA Tech_3 AGHS4_LA Tech_4

Number of Observations

Number of Observations Read 144
Number of Observations Used 144

Dimensions

Dimensions
G-side Cov. Parameters 2
R-side Cov. Parameters 1
Columns in X 84
Columns in Z per Subject 1
Subjects (Blocks in V) 48
Max Obs per Subject 3

Optimization Information

Optimization Information
Optimization Technique Dual Quasi-Newton
Parameters in Optimization 2
Lower Boundaries 2
Upper Boundaries 1
Fixed Effects Profiled
Residual Variance Profiled
Starting From Data

Iteration History

Iteration History
Iteration Restarts Evaluations Objective
Function
Change Max
Gradient
0 0 4 640.72210261 . 1.46E-13

Convergence Status

Convergence criterion (ABSGCONV=0.00001) satisfied.

Hessian Matrix

Hessian Matrix
Cov Parm Subject Row Col1 Col2
Variance GLB 1 80.7379  
AR(1) GLB 2    

Fit Statistics

Fit Statistics
-2 Res Log Likelihood 640.72
AIC (smaller is better) 646.72
AICC (smaller is better) 646.95
BIC (smaller is better) 652.34
CAIC (smaller is better) 655.34
HQIC (smaller is better) 648.84
Generalized Chi-Square 1275.07
Gener. Chi-Square / DF 11.81

Covariance Parameter Estimates

Covariance Parameter Estimates
Cov Parm Subject Estimate Standard
Error
Variance GLB 2.5015 1.6529
AR(1) GLB 0 .
Residual   11.8062 1.9677

Solutions for Fixed Effects

Solutions for Fixed Effects
Effect Year Location Family Estimate Standard
Error
DF t Value Pr > |t|
Intercept       48.3072 1.8913 101.8 25.54 <.0001
Year First     -34.9458 2.4296 72 -14.38 <.0001
Year Initial     -43.5014 2.4296 72 -17.90 <.0001
Year Second     0 . . . .
Location   Hill Farm   -8.7071 2.6747 101.8 -3.26 0.0015
Location   LA Tech   0 . . . .
Year*Location First Hill Farm   9.0771 3.4360 72 2.64 0.0101
Year*Location First LA Tech   0 . . . .
Year*Location Initial Hill Farm   9.6024 3.4360 72 2.79 0.0067
Year*Location Initial LA Tech   0 . . . .
Year*Location Second Hill Farm   0 . . . .
Year*Location Second LA Tech   0 . . . .
Family     AGH2 -12.1855 2.6747 101.8 -4.56 <.0001
Family     AGH25 -6.4656 2.6747 101.8 -2.42 0.0174
Family     AGHS1 -0.8208 2.6747 101.8 -0.31 0.7596
Family     AGHS2 2.9486 2.6747 101.8 1.10 0.2729
Family     AGHS3 -1.5495 2.6747 101.8 -0.58 0.5637
Family     AGHS4 0 . . . .
Year*Family First   AGH2 9.1584 3.4360 72 2.67 0.0095
Year*Family First   AGH25 3.8052 3.4360 72 1.11 0.2718
Year*Family First   AGHS1 -0.1953 3.4360 72 -0.06 0.9548
Year*Family First   AGHS2 -2.9219 3.4360 72 -0.85 0.3979
Year*Family First   AGHS3 0.4433 3.4360 72 0.13 0.8977
Year*Family First   AGHS4 0 . . . .
Year*Family Initial   AGH2 12.1499 3.4360 72 3.54 0.0007
Year*Family Initial   AGH25 6.5162 3.4360 72 1.90 0.0619
Year*Family Initial   AGHS1 0.5714 3.4360 72 0.17 0.8684
Year*Family Initial   AGHS2 -3.2294 3.4360 72 -0.94 0.3504
Year*Family Initial   AGHS3 1.6059 3.4360 72 0.47 0.6417
Year*Family Initial   AGHS4 0 . . . .
Year*Family Second   AGH2 0 . . . .
Year*Family Second   AGH25 0 . . . .
Year*Family Second   AGHS1 0 . . . .
Year*Family Second   AGHS2 0 . . . .
Year*Family Second   AGHS3 0 . . . .
Year*Family Second   AGHS4 0 . . . .
Location*Family   Hill Farm AGH2 6.2013 3.7826 101.8 1.64 0.1042
Location*Family   Hill Farm AGH25 1.1020 3.7826 101.8 0.29 0.7714
Location*Family   Hill Farm AGHS1 2.9179 3.7826 101.8 0.77 0.4423
Location*Family   Hill Farm AGHS2 -3.6520 3.7826 101.8 -0.97 0.3366
Location*Family   Hill Farm AGHS3 -1.4682 3.7826 101.8 -0.39 0.6987
Location*Family   Hill Farm AGHS4 0 . . . .
Location*Family   LA Tech AGH2 0 . . . .
Location*Family   LA Tech AGH25 0 . . . .
Location*Family   LA Tech AGHS1 0 . . . .
Location*Family   LA Tech AGHS2 0 . . . .
Location*Family   LA Tech AGHS3 0 . . . .
Location*Family   LA Tech AGHS4 0 . . . .
Year*Location*Family First Hill Farm AGH2 -4.4378 4.8593 72 -0.91 0.3641
Year*Location*Family First Hill Farm AGH25 1.3168 4.8593 72 0.27 0.7872
Year*Location*Family First Hill Farm AGHS1 0.3735 4.8593 72 0.08 0.9389
Year*Location*Family First Hill Farm AGHS2 5.2515 4.8593 72 1.08 0.2834
Year*Location*Family First Hill Farm AGHS3 3.3401 4.8593 72 0.69 0.4941
Year*Location*Family First Hill Farm AGHS4 0 . . . .
Year*Location*Family First LA Tech AGH2 0 . . . .
Year*Location*Family First LA Tech AGH25 0 . . . .
Year*Location*Family First LA Tech AGHS1 0 . . . .
Year*Location*Family First LA Tech AGHS2 0 . . . .
Year*Location*Family First LA Tech AGHS3 0 . . . .
Year*Location*Family First LA Tech AGHS4 0 . . . .
Year*Location*Family Initial Hill Farm AGH2 -6.4171 4.8593 72 -1.32 0.1908
Year*Location*Family Initial Hill Farm AGH25 -1.7337 4.8593 72 -0.36 0.7223
Year*Location*Family Initial Hill Farm AGHS1 -2.6938 4.8593 72 -0.55 0.5810
Year*Location*Family Initial Hill Farm AGHS2 3.1032 4.8593 72 0.64 0.5251
Year*Location*Family Initial Hill Farm AGHS3 1.0562 4.8593 72 0.22 0.8285
Year*Location*Family Initial Hill Farm AGHS4 0 . . . .
Year*Location*Family Initial LA Tech AGH2 0 . . . .
Year*Location*Family Initial LA Tech AGH25 0 . . . .
Year*Location*Family Initial LA Tech AGHS1 0 . . . .
Year*Location*Family Initial LA Tech AGHS2 0 . . . .
Year*Location*Family Initial LA Tech AGHS3 0 . . . .
Year*Location*Family Initial LA Tech AGHS4 0 . . . .
Year*Location*Family Second Hill Farm AGH2 0 . . . .
Year*Location*Family Second Hill Farm AGH25 0 . . . .
Year*Location*Family Second Hill Farm AGHS1 0 . . . .
Year*Location*Family Second Hill Farm AGHS2 0 . . . .
Year*Location*Family Second Hill Farm AGHS3 0 . . . .
Year*Location*Family Second Hill Farm AGHS4 0 . . . .
Year*Location*Family Second LA Tech AGH2 0 . . . .
Year*Location*Family Second LA Tech AGH25 0 . . . .
Year*Location*Family Second LA Tech AGHS1 0 . . . .
Year*Location*Family Second LA Tech AGHS2 0 . . . .
Year*Location*Family Second LA Tech AGHS3 0 . . . .
Year*Location*Family Second LA Tech AGHS4 0 . . . .

Type III Tests of Fixed Effects

Type III Tests of Fixed Effects
Effect Num DF Den DF F Value Pr > F
Year 2 72 1477.65 <.0001
Location 1 36 5.24 0.0280
Year*Location 2 72 29.73 <.0001
Family 5 36 3.78 0.0075
Year*Family 10 72 3.10 0.0025
Location*Family 5 36 0.56 0.7268
Year*Location*Family 10 72 0.65 0.7634

Solution for Random Effects

Solution for Random Effects
Effect Subject Estimate Std Err Pred DF t Value Pr > |t| Alpha Lower Upper
Intercept GLB AGH25_Hill Farm_1 1.0602 1.4398 8.699 0.74 0.4809 0.05 -2.2142 4.3346
Intercept GLB AGH25_Hill Farm_2 -0.05423 1.4398 8.699 -0.04 0.9708 0.05 -3.3287 3.2202
Intercept GLB AGH25_Hill Farm_3 -1.0796 1.4398 8.699 -0.75 0.4732 0.05 -4.3540 2.1949
Intercept GLB AGH25_Hill Farm_4 0.07359 1.4398 8.699 0.05 0.9604 0.05 -3.2008 3.3480
Intercept GLB AGH25_LA Tech_1 -0.5470 1.4398 8.699 -0.38 0.7131 0.05 -3.8215 2.7274
Intercept GLB AGH25_LA Tech_2 0.09153 1.4398 8.699 0.06 0.9508 0.05 -3.1829 3.3660
Intercept GLB AGH25_LA Tech_3 -0.6572 1.4398 8.699 -0.46 0.6593 0.05 -3.9316 2.6172
Intercept GLB AGH25_LA Tech_4 1.1127 1.4398 8.699 0.77 0.4601 0.05 -2.1617 4.3872
Intercept GLB AGH2_Hill Farm_1 0.3198 1.4398 8.699 0.22 0.8293 0.05 -2.9546 3.5943
Intercept GLB AGH2_Hill Farm_2 0.3363 1.4398 8.699 0.23 0.8207 0.05 -2.9381 3.6108
Intercept GLB AGH2_Hill Farm_3 0.1483 1.4398 8.699 0.10 0.9203 0.05 -3.1261 3.4228
Intercept GLB AGH2_Hill Farm_4 -0.8045 1.4398 8.699 -0.56 0.5904 0.05 -4.0789 2.4699
Intercept GLB AGH2_LA Tech_1 0.2711 1.4398 8.699 0.19 0.8550 0.05 -3.0033 3.5456
Intercept GLB AGH2_LA Tech_2 -0.1484 1.4398 8.699 -0.10 0.9203 0.05 -3.4228 3.1261
Intercept GLB AGH2_LA Tech_3 0.6397 1.4398 8.699 0.44 0.6677 0.05 -2.6348 3.9141
Intercept GLB AGH2_LA Tech_4 -0.7625 1.4398 8.699 -0.53 0.6097 0.05 -4.0369 2.5120
Intercept GLB AGHS1_Hill Farm_1 0.7940 1.4398 8.699 0.55 0.5952 0.05 -2.4805 4.0684
Intercept GLB AGHS1_Hill Farm_2 -0.4143 1.4398 8.699 -0.29 0.7803 0.05 -3.6888 2.8601
Intercept GLB AGHS1_Hill Farm_3 0.7606 1.4398 8.699 0.53 0.6105 0.05 -2.5138 4.0350
Intercept GLB AGHS1_Hill Farm_4 -1.1402 1.4398 8.699 -0.79 0.4495 0.05 -4.4146 2.1342
Intercept GLB AGHS1_LA Tech_1 -0.6798 1.4398 8.699 -0.47 0.6484 0.05 -3.9543 2.5946
Intercept GLB AGHS1_LA Tech_2 1.3655 1.4398 8.699 0.95 0.3686 0.05 -1.9090 4.6399
Intercept GLB AGHS1_LA Tech_3 -0.1045 1.4398 8.699 -0.07 0.9438 0.05 -3.3790 3.1699
Intercept GLB AGHS1_LA Tech_4 -0.5811 1.4398 8.699 -0.40 0.6963 0.05 -3.8555 2.6933
Intercept GLB AGHS2_Hill Farm_1 -1.2741 1.4398 8.699 -0.88 0.4000 0.05 -4.5485 2.0003
Intercept GLB AGHS2_Hill Farm_2 -1.9036 1.4398 8.699 -1.32 0.2199 0.05 -5.1780 1.3709
Intercept GLB AGHS2_Hill Farm_3 2.0137 1.4398 8.699 1.40 0.1966 0.05 -1.2607 5.2881
Intercept GLB AGHS2_Hill Farm_4 1.1639 1.4398 8.699 0.81 0.4404 0.05 -2.1105 4.4384
Intercept GLB AGHS2_LA Tech_1 -1.1906 1.4398 8.699 -0.83 0.4304 0.05 -4.4650 2.0838
Intercept GLB AGHS2_LA Tech_2 -0.4809 1.4398 8.699 -0.33 0.7463 0.05 -3.7553 2.7935
Intercept GLB AGHS2_LA Tech_3 0.5704 1.4398 8.699 0.40 0.7015 0.05 -2.7040 3.8448
Intercept GLB AGHS2_LA Tech_4 1.1011 1.4398 8.699 0.76 0.4647 0.05 -2.1733 4.3755
Intercept GLB AGHS3_Hill Farm_1 1.5679 1.4398 8.699 1.09 0.3054 0.05 -1.7066 4.8423
Intercept GLB AGHS3_Hill Farm_2 -1.5387 1.4398 8.699 -1.07 0.3140 0.05 -4.8132 1.7357
Intercept GLB AGHS3_Hill Farm_3 0.3698 1.4398 8.699 0.26 0.8033 0.05 -2.9046 3.6442
Intercept GLB AGHS3_Hill Farm_4 -0.3989 1.4398 8.699 -0.28 0.7882 0.05 -3.6734 2.8755
Intercept GLB AGHS3_LA Tech_1 -0.01348 1.4398 8.699 -0.01 0.9927 0.05 -3.2879 3.2609
Intercept GLB AGHS3_LA Tech_2 -0.1444 1.4398 8.699 -0.10 0.9224 0.05 -3.4188 3.1301
Intercept GLB AGHS3_LA Tech_3 -0.3383 1.4398 8.699 -0.23 0.8197 0.05 -3.6127 2.9361
Intercept GLB AGHS3_LA Tech_4 0.4961 1.4398 8.699 0.34 0.7386 0.05 -2.7783 3.7706
Intercept GLB AGHS4_Hill Farm_1 -0.4645 1.4398 8.699 -0.32 0.7546 0.05 -3.7389 2.8100
Intercept GLB AGHS4_Hill Farm_2 1.3958 1.4398 8.699 0.97 0.3585 0.05 -1.8786 4.6702
Intercept GLB AGHS4_Hill Farm_3 -1.1545 1.4398 8.699 -0.80 0.4440 0.05 -4.4289 2.1199
Intercept GLB AGHS4_Hill Farm_4 0.2232 1.4398 8.699 0.16 0.8804 0.05 -3.0512 3.4976
Intercept GLB AGHS4_LA Tech_1 -0.5897 1.4398 8.699 -0.41 0.6920 0.05 -3.8641 2.6847
Intercept GLB AGHS4_LA Tech_2 0.08132 1.4398 8.699 0.06 0.9562 0.05 -3.1931 3.3558
Intercept GLB AGHS4_LA Tech_3 0.1598 1.4398 8.699 0.11 0.9142 0.05 -3.1147 3.4342
Intercept GLB AGHS4_LA Tech_4 0.3486 1.4398 8.699 0.24 0.8143 0.05 -2.9258 3.6230

Slices

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year First 1 101.8 4.04 0.0471

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Location _Location Estimate Standard Error DF t Value Pr > |t| Adj P
Year First Hill Farm LA Tech 2.1942 1.0919 101.8 2.01 0.0471 0.0482

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Location Estimate  
Year First Hill Farm 14.2584 A
Year First      
Year First LA Tech 12.0642 B

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Initial 1 101.8 0.33 0.5645

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Location _Location Estimate Standard Error DF t Value Pr > |t| Adj P
Year Initial Hill Farm LA Tech 0.6313 1.0919 101.8 0.58 0.5645 0.5650

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Location Estimate  
Year Initial Hill Farm 5.3606 A
Year Initial     A
Year Initial LA Tech 4.7294 A

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Second 1 101.8 51.77 <.0001

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Location _Location Estimate Standard Error DF t Value Pr > |t| Adj P
Year Second Hill Farm LA Tech -7.8569 1.0919 101.8 -7.20 <.0001 <.0001

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Location Estimate  
Year Second LA Tech 45.2951 A
Year Second      
Year Second Hill Farm 37.4382 B

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year First 5 101.8 0.77 0.5702

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Family _Family Estimate Standard Error DF t Value Pr > |t| Adj P
Year First AGH2 AGH25 -0.6944 1.8913 101.8 -0.37 0.7143 0.9991
Year First AGH2 AGHS1 -2.7749 1.8913 101.8 -1.47 0.1454 0.6859
Year First AGH2 AGHS2 -2.9718 1.8913 101.8 -1.57 0.1192 0.6198
Year First AGH2 AGHS3 -1.9752 1.8913 101.8 -1.04 0.2988 0.9012
Year First AGH2 AGHS4 -2.1454 1.8913 101.8 -1.13 0.2593 0.8654
Year First AGH25 AGHS1 -2.0806 1.8913 101.8 -1.10 0.2739 0.8798
Year First AGH25 AGHS2 -2.2774 1.8913 101.8 -1.20 0.2313 0.8334
Year First AGH25 AGHS3 -1.2808 1.8913 101.8 -0.68 0.4998 0.9839
Year First AGH25 AGHS4 -1.4510 1.8913 101.8 -0.77 0.4447 0.9721
Year First AGHS1 AGHS2 -0.1969 1.8913 101.8 -0.10 0.9173 1.0000
Year First AGHS1 AGHS3 0.7998 1.8913 101.8 0.42 0.6733 0.9982
Year First AGHS1 AGHS4 0.6296 1.8913 101.8 0.33 0.7399 0.9994
Year First AGHS2 AGHS3 0.9966 1.8913 101.8 0.53 0.5994 0.9949
Year First AGHS2 AGHS4 0.8264 1.8913 101.8 0.44 0.6631 0.9979
Year First AGHS3 AGHS4 -0.1702 1.8913 101.8 -0.09 0.9285 1.0000

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Family Estimate  
Year First AGHS2 14.3728 A
Year First     A
Year First AGHS1 14.1760 A
Year First     A
Year First AGHS4 13.5464 A
Year First     A
Year First AGHS3 13.3762 A
Year First     A
Year First AGH25 12.0954 A
Year First     A
Year First AGH2 11.4010 A

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Initial 5 101.8 0.02 0.9998

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Family _Family Estimate Standard Error DF t Value Pr > |t| Adj P
Year Initial AGH2 AGH25 0.1218 1.8913 101.8 0.06 0.9488 1.0000
Year Initial AGH2 AGHS1 -0.00609 1.8913 101.8 -0.00 0.9974 1.0000
Year Initial AGH2 AGHS2 0.4117 1.8913 101.8 0.22 0.8281 0.9999
Year Initial AGH2 AGHS3 0.006111 1.8913 101.8 0.00 0.9974 1.0000
Year Initial AGH2 AGHS4 -0.1435 1.8913 101.8 -0.08 0.9397 1.0000
Year Initial AGH25 AGHS1 -0.1279 1.8913 101.8 -0.07 0.9462 1.0000
Year Initial AGH25 AGHS2 0.2899 1.8913 101.8 0.15 0.8785 1.0000
Year Initial AGH25 AGHS3 -0.1157 1.8913 101.8 -0.06 0.9513 1.0000
Year Initial AGH25 AGHS4 -0.2653 1.8913 101.8 -0.14 0.8887 1.0000
Year Initial AGHS1 AGHS2 0.4178 1.8913 101.8 0.22 0.8256 0.9999
Year Initial AGHS1 AGHS3 0.01220 1.8913 101.8 0.01 0.9949 1.0000
Year Initial AGHS1 AGHS4 -0.1374 1.8913 101.8 -0.07 0.9422 1.0000
Year Initial AGHS2 AGHS3 -0.4056 1.8913 101.8 -0.21 0.8306 0.9999
Year Initial AGHS2 AGHS4 -0.5552 1.8913 101.8 -0.29 0.7697 0.9997
Year Initial AGHS3 AGHS4 -0.1496 1.8913 101.8 -0.08 0.9371 1.0000

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Family Estimate  
Year Initial AGHS4 5.2535 A
Year Initial     A
Year Initial AGHS1 5.1161 A
Year Initial     A
Year Initial AGH2 5.1100 A
Year Initial     A
Year Initial AGHS3 5.1039 A
Year Initial     A
Year Initial AGH25 4.9882 A
Year Initial     A
Year Initial AGHS2 4.6983 A

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Second 5 101.8 9.42 <.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Family _Family Estimate Standard Error DF t Value Pr > |t| Adj P
Year Second AGH2 AGH25 -3.1702 1.8913 101.8 -1.68 0.0968 0.5517
Year Second AGH2 AGHS1 -9.7229 1.8913 101.8 -5.14 <.0001 <.0001
Year Second AGH2 AGHS2 -10.2074 1.8913 101.8 -5.40 <.0001 <.0001
Year Second AGH2 AGHS3 -6.8012 1.8913 101.8 -3.60 0.0005 0.0075
Year Second AGH2 AGHS4 -9.0848 1.8913 101.8 -4.80 <.0001 0.0001
Year Second AGH25 AGHS1 -6.5527 1.8913 101.8 -3.46 0.0008 0.0112
Year Second AGH25 AGHS2 -7.0372 1.8913 101.8 -3.72 0.0003 0.0051
Year Second AGH25 AGHS3 -3.6310 1.8913 101.8 -1.92 0.0577 0.3987
Year Second AGH25 AGHS4 -5.9146 1.8913 101.8 -3.13 0.0023 0.0294
Year Second AGHS1 AGHS2 -0.4845 1.8913 101.8 -0.26 0.7983 0.9998
Year Second AGHS1 AGHS3 2.9217 1.8913 101.8 1.54 0.1255 0.6368
Year Second AGHS1 AGHS4 0.6381 1.8913 101.8 0.34 0.7365 0.9994
Year Second AGHS2 AGHS3 3.4062 1.8913 101.8 1.80 0.0747 0.4716
Year Second AGHS2 AGHS4 1.1226 1.8913 101.8 0.59 0.5541 0.9912
Year Second AGHS3 AGHS4 -2.2836 1.8913 101.8 -1.21 0.2301 0.8319

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Family Estimate  
Year Second AGHS2 45.0763   A
Year Second       A
Year Second AGHS1 44.5918   A
Year Second       A
Year Second AGHS4 43.9537   A
Year Second       A
Year Second AGHS3 41.6701 B A
Year Second     B  
Year Second AGH25 38.0391 B C
Year Second       C
Year Second AGH2 34.8689   C

LS-Means Results

Family LS-Means

Family Least Squares Means
Family Estimate Standard
Error
DF t Value Pr > |t|
AGH2 17.1266 0.8970 36 19.09 <.0001
AGH25 18.3742 0.8970 36 20.48 <.0001
AGHS1 21.2946 0.8970 36 23.74 <.0001
AGHS2 21.3825 0.8970 36 23.84 <.0001
AGHS3 20.0501 0.8970 36 22.35 <.0001
AGHS4 20.9179 0.8970 36 23.32 <.0001

Family Diffs

Differences of Family Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Family _Family Estimate Standard Error DF t Value Pr > |t| Adj P
AGH2 AGH25 -1.2476 1.2686 36 -0.98 0.3319 0.9202
AGH2 AGHS1 -4.1680 1.2686 36 -3.29 0.0023 0.0255
AGH2 AGHS2 -4.2558 1.2686 36 -3.35 0.0019 0.0214
AGH2 AGHS3 -2.9234 1.2686 36 -2.30 0.0271 0.2187
AGH2 AGHS4 -3.7912 1.2686 36 -2.99 0.0050 0.0524
AGH25 AGHS1 -2.9204 1.2686 36 -2.30 0.0272 0.2196
AGH25 AGHS2 -3.0082 1.2686 36 -2.37 0.0232 0.1934
AGH25 AGHS3 -1.6758 1.2686 36 -1.32 0.1948 0.7716
AGH25 AGHS4 -2.5436 1.2686 36 -2.01 0.0525 0.3594
AGHS1 AGHS2 -0.08785 1.2686 36 -0.07 0.9452 1.0000
AGHS1 AGHS3 1.2445 1.2686 36 0.98 0.3331 0.9210
AGHS1 AGHS4 0.3768 1.2686 36 0.30 0.7682 0.9997
AGHS2 AGHS3 1.3324 1.2686 36 1.05 0.3006 0.8974
AGHS2 AGHS4 0.4646 1.2686 36 0.37 0.7163 0.9991
AGHS3 AGHS4 -0.8678 1.2686 36 -0.68 0.4983 0.9826

Family Lines

Tukey Grouping for Family Least Squares Means (Alpha=0.05)
LS-means with the same letter are not significantly different.
Family Estimate  
AGHS2 21.3825   A
      A
AGHS1 21.2946   A
      A
AGHS4 20.9179 B A
    B A
AGHS3 20.0501 B A
    B A
AGH25 18.3742 B A
    B  
AGH2 17.1266 B  

SAS Output DBH

htmltools::includeHTML("https://raw.githubusercontent.com/azd169/SAS_outputs/main/hane-results-DBH.html")
## Warning: `includeHTML()` was provided a `path` that appears to be a complete HTML document.
## ✖ Path: https://raw.githubusercontent.com/azd169/SAS_outputs/main/hane-results-DBH.html
## ℹ Use `tags$iframe()` to include an HTML document. You can either ensure `path` is accessible in your app or document (see e.g. `shiny::addResourcePath()`) and pass the relative path to the `src` argument. Or you can read the contents of `path` and pass the contents to `srcdoc`.
Results: hane.sas

Results: hane.sas

The GLIMMIX Procedure

The Glimmix Procedure

Model Information

Model Information
Data Set WORK.DBH
Response Variable DBH
Response Distribution Gaussian
Link Function Identity
Variance Function Default
Variance Matrix Blocked By GLB
Estimation Technique Restricted Maximum Likelihood
Degrees of Freedom Method Kenward-Roger
Fixed Effects SE Adjustment Kenward-Roger

Class Level Information

Class Level Information
Class Levels Values
Year 5 Fifth Fourth Seventh Sixth Third
Location 2 Hill Farm LA Tech
Block 4 1 2 3 4
Family 6 AGH2 AGH25 AGHS1 AGHS2 AGHS3 AGHS4
GLB 48 AGH25_Hill Farm_1 AGH25_Hill Farm_2 AGH25_Hill Farm_3 AGH25_Hill Farm_4 AGH25_LA Tech_1 AGH25_LA Tech_2 AGH25_LA Tech_3 AGH25_LA Tech_4 AGH2_Hill Farm_1 AGH2_Hill Farm_2 AGH2_Hill Farm_3 AGH2_Hill Farm_4 AGH2_LA Tech_1 AGH2_LA Tech_2 AGH2_LA Tech_3 AGH2_LA Tech_4 AGHS1_Hill Farm_1 AGHS1_Hill Farm_2 AGHS1_Hill Farm_3 AGHS1_Hill Farm_4 AGHS1_LA Tech_1 AGHS1_LA Tech_2 AGHS1_LA Tech_3 AGHS1_LA Tech_4 AGHS2_Hill Farm_1 AGHS2_Hill Farm_2 AGHS2_Hill Farm_3 AGHS2_Hill Farm_4 AGHS2_LA Tech_1 AGHS2_LA Tech_2 AGHS2_LA Tech_3 AGHS2_LA Tech_4 AGHS3_Hill Farm_1 AGHS3_Hill Farm_2 AGHS3_Hill Farm_3 AGHS3_Hill Farm_4 AGHS3_LA Tech_1 AGHS3_LA Tech_2 AGHS3_LA Tech_3 AGHS3_LA Tech_4 AGHS4_Hill Farm_1 AGHS4_Hill Farm_2 AGHS4_Hill Farm_3 AGHS4_Hill Farm_4 AGHS4_LA Tech_1 AGHS4_LA Tech_2 AGHS4_LA Tech_3 AGHS4_LA Tech_4

Number of Observations

Number of Observations Read 240
Number of Observations Used 240

Dimensions

Dimensions
G-side Cov. Parameters 2
R-side Cov. Parameters 1
Columns in X 126
Columns in Z per Subject 1
Subjects (Blocks in V) 48
Max Obs per Subject 5

Optimization Information

Optimization Information
Optimization Technique Dual Quasi-Newton
Parameters in Optimization 2
Lower Boundaries 2
Upper Boundaries 1
Fixed Effects Profiled
Residual Variance Profiled
Starting From Data

Iteration History

Iteration History
Iteration Restarts Evaluations Objective
Function
Change Max
Gradient
0 0 4 337.79560627 . 5.28E-13

Convergence Status

Convergence criterion (ABSGCONV=0.00001) satisfied.

Hessian Matrix

Hessian Matrix
Cov Parm Subject Row Col1 Col2
Variance GLB 1 7.5640  
AR(1) GLB 2    

Fit Statistics

Fit Statistics
-2 Res Log Likelihood 337.80
AIC (smaller is better) 343.80
AICC (smaller is better) 343.93
BIC (smaller is better) 349.41
CAIC (smaller is better) 352.41
HQIC (smaller is better) 345.92
Generalized Chi-Square 27.50
Gener. Chi-Square / DF 0.15

Covariance Parameter Estimates

Covariance Parameter Estimates
Cov Parm Subject Estimate Standard
Error
Variance GLB 0.2675 0.07035
AR(1) GLB 0 .
Residual   0.1528 0.01800

Solutions for Fixed Effects

Solutions for Fixed Effects
Effect Year Location Family Estimate Standard
Error
DF t Value Pr > |t|
Intercept       3.8518 0.3241 68.68 11.88 <.0001
Year Fifth     6.8444 0.2764 144 24.77 <.0001
Year Fourth     3.7379 0.2764 144 13.53 <.0001
Year Seventh     9.9398 0.2764 144 35.97 <.0001
Year Sixth     8.7141 0.2764 144 31.53 <.0001
Year Third     0 . . . .
Location   Hill Farm   -0.9135 0.4584 68.68 -1.99 0.0502
Location   LA Tech   0 . . . .
Year*Location Fifth Hill Farm   -1.2948 0.3908 144 -3.31 0.0012
Year*Location Fifth LA Tech   0 . . . .
Year*Location Fourth Hill Farm   -0.5688 0.3908 144 -1.46 0.1477
Year*Location Fourth LA Tech   0 . . . .
Year*Location Seventh Hill Farm   -0.6350 0.3908 144 -1.62 0.1064
Year*Location Seventh LA Tech   0 . . . .
Year*Location Sixth Hill Farm   -0.4200 0.3908 144 -1.07 0.2844
Year*Location Sixth LA Tech   0 . . . .
Year*Location Third Hill Farm   0 . . . .
Year*Location Third LA Tech   0 . . . .
Family     AGH2 -2.1578 0.4584 68.68 -4.71 <.0001
Family     AGH25 -1.9694 0.4584 68.68 -4.30 <.0001
Family     AGHS1 -0.3172 0.4584 68.68 -0.69 0.4913
Family     AGHS2 0.1063 0.4584 68.68 0.23 0.8173
Family     AGHS3 -1.2254 0.4584 68.68 -2.67 0.0094
Family     AGHS4 0 . . . .
Year*Family Fifth   AGH2 -0.5357 0.3908 144 -1.37 0.1726
Year*Family Fifth   AGH25 -0.2856 0.3908 144 -0.73 0.4661
Year*Family Fifth   AGHS1 -1.4184 0.3908 144 -3.63 0.0004
Year*Family Fifth   AGHS2 -1.6092 0.3908 144 -4.12 <.0001
Year*Family Fifth   AGHS3 -0.1346 0.3908 144 -0.34 0.7311
Year*Family Fifth   AGHS4 0 . . . .
Year*Family Fourth   AGH2 -0.3221 0.3908 144 -0.82 0.4112
Year*Family Fourth   AGH25 -0.01489 0.3908 144 -0.04 0.9697
Year*Family Fourth   AGHS1 -0.6544 0.3908 144 -1.67 0.0962
Year*Family Fourth   AGHS2 -0.7991 0.3908 144 -2.04 0.0427
Year*Family Fourth   AGHS3 0.06217 0.3908 144 0.16 0.8738
Year*Family Fourth   AGHS4 0 . . . .
Year*Family Seventh   AGH2 -0.04498 0.3908 144 -0.12 0.9085
Year*Family Seventh   AGH25 -0.2542 0.3908 144 -0.65 0.5165
Year*Family Seventh   AGHS1 -1.5222 0.3908 144 -3.89 0.0002
Year*Family Seventh   AGHS2 -1.6619 0.3908 144 -4.25 <.0001
Year*Family Seventh   AGHS3 -0.1730 0.3908 144 -0.44 0.6586
Year*Family Seventh   AGHS4 0 . . . .
Year*Family Sixth   AGH2 -0.4924 0.3908 144 -1.26 0.2097
Year*Family Sixth   AGH25 -0.2868 0.3908 144 -0.73 0.4643
Year*Family Sixth   AGHS1 -1.2209 0.3908 144 -3.12 0.0022
Year*Family Sixth   AGHS2 -1.6868 0.3908 144 -4.32 <.0001
Year*Family Sixth   AGHS3 0.09563 0.3908 144 0.24 0.8070
Year*Family Sixth   AGHS4 0 . . . .
Year*Family Third   AGH2 0 . . . .
Year*Family Third   AGH25 0 . . . .
Year*Family Third   AGHS1 0 . . . .
Year*Family Third   AGHS2 0 . . . .
Year*Family Third   AGHS3 0 . . . .
Year*Family Third   AGHS4 0 . . . .
Location*Family   Hill Farm AGH2 0.7712 0.6483 68.68 1.19 0.2383
Location*Family   Hill Farm AGH25 0.5841 0.6483 68.68 0.90 0.3708
Location*Family   Hill Farm AGHS1 0.2693 0.6483 68.68 0.42 0.6792
Location*Family   Hill Farm AGHS2 -0.4456 0.6483 68.68 -0.69 0.4941
Location*Family   Hill Farm AGHS3 0.2826 0.6483 68.68 0.44 0.6642
Location*Family   Hill Farm AGHS4 0 . . . .
Location*Family   LA Tech AGH2 0 . . . .
Location*Family   LA Tech AGH25 0 . . . .
Location*Family   LA Tech AGHS1 0 . . . .
Location*Family   LA Tech AGHS2 0 . . . .
Location*Family   LA Tech AGHS3 0 . . . .
Location*Family   LA Tech AGHS4 0 . . . .
Year*Location*Family Fifth Hill Farm AGH2 -0.4343 0.5527 144 -0.79 0.4333
Year*Location*Family Fifth Hill Farm AGH25 -0.2834 0.5527 144 -0.51 0.6089
Year*Location*Family Fifth Hill Farm AGHS1 0.06586 0.5527 144 0.12 0.9053
Year*Location*Family Fifth Hill Farm AGHS2 0.7027 0.5527 144 1.27 0.2057
Year*Location*Family Fifth Hill Farm AGHS3 -0.1706 0.5527 144 -0.31 0.7580
Year*Location*Family Fifth Hill Farm AGHS4 0 . . . .
Year*Location*Family Fifth LA Tech AGH2 0 . . . .
Year*Location*Family Fifth LA Tech AGH25 0 . . . .
Year*Location*Family Fifth LA Tech AGHS1 0 . . . .
Year*Location*Family Fifth LA Tech AGHS2 0 . . . .
Year*Location*Family Fifth LA Tech AGHS3 0 . . . .
Year*Location*Family Fifth LA Tech AGHS4 0 . . . .
Year*Location*Family Fourth Hill Farm AGH2 -0.9294 0.5527 144 -1.68 0.0948
Year*Location*Family Fourth Hill Farm AGH25 -1.0055 0.5527 144 -1.82 0.0710
Year*Location*Family Fourth Hill Farm AGHS1 -0.2051 0.5527 144 -0.37 0.7112
Year*Location*Family Fourth Hill Farm AGHS2 -0.04451 0.5527 144 -0.08 0.9359
Year*Location*Family Fourth Hill Farm AGHS3 -0.7046 0.5527 144 -1.27 0.2044
Year*Location*Family Fourth Hill Farm AGHS4 0 . . . .
Year*Location*Family Fourth LA Tech AGH2 0 . . . .
Year*Location*Family Fourth LA Tech AGH25 0 . . . .
Year*Location*Family Fourth LA Tech AGHS1 0 . . . .
Year*Location*Family Fourth LA Tech AGHS2 0 . . . .
Year*Location*Family Fourth LA Tech AGHS3 0 . . . .
Year*Location*Family Fourth LA Tech AGHS4 0 . . . .
Year*Location*Family Seventh Hill Farm AGH2 -0.8942 0.5527 144 -1.62 0.1079
Year*Location*Family Seventh Hill Farm AGH25 -0.6181 0.5527 144 -1.12 0.2653
Year*Location*Family Seventh Hill Farm AGHS1 -0.3015 0.5527 144 -0.55 0.5862
Year*Location*Family Seventh Hill Farm AGHS2 0.4281 0.5527 144 0.77 0.4398
Year*Location*Family Seventh Hill Farm AGHS3 -0.8381 0.5527 144 -1.52 0.1316
Year*Location*Family Seventh Hill Farm AGHS4 0 . . . .
Year*Location*Family Seventh LA Tech AGH2 0 . . . .
Year*Location*Family Seventh LA Tech AGH25 0 . . . .
Year*Location*Family Seventh LA Tech AGHS1 0 . . . .
Year*Location*Family Seventh LA Tech AGHS2 0 . . . .
Year*Location*Family Seventh LA Tech AGHS3 0 . . . .
Year*Location*Family Seventh LA Tech AGHS4 0 . . . .
Year*Location*Family Sixth Hill Farm AGH2 -0.6702 0.5527 144 -1.21 0.2273
Year*Location*Family Sixth Hill Farm AGH25 -0.7039 0.5527 144 -1.27 0.2049
Year*Location*Family Sixth Hill Farm AGHS1 -0.5868 0.5527 144 -1.06 0.2902
Year*Location*Family Sixth Hill Farm AGHS2 0.3045 0.5527 144 0.55 0.5825
Year*Location*Family Sixth Hill Farm AGHS3 -0.5005 0.5527 144 -0.91 0.3667
Year*Location*Family Sixth Hill Farm AGHS4 0 . . . .
Year*Location*Family Sixth LA Tech AGH2 0 . . . .
Year*Location*Family Sixth LA Tech AGH25 0 . . . .
Year*Location*Family Sixth LA Tech AGHS1 0 . . . .
Year*Location*Family Sixth LA Tech AGHS2 0 . . . .
Year*Location*Family Sixth LA Tech AGHS3 0 . . . .
Year*Location*Family Sixth LA Tech AGHS4 0 . . . .
Year*Location*Family Third Hill Farm AGH2 0 . . . .
Year*Location*Family Third Hill Farm AGH25 0 . . . .
Year*Location*Family Third Hill Farm AGHS1 0 . . . .
Year*Location*Family Third Hill Farm AGHS2 0 . . . .
Year*Location*Family Third Hill Farm AGHS3 0 . . . .
Year*Location*Family Third Hill Farm AGHS4 0 . . . .
Year*Location*Family Third LA Tech AGH2 0 . . . .
Year*Location*Family Third LA Tech AGH25 0 . . . .
Year*Location*Family Third LA Tech AGHS1 0 . . . .
Year*Location*Family Third LA Tech AGHS2 0 . . . .
Year*Location*Family Third LA Tech AGHS3 0 . . . .
Year*Location*Family Third LA Tech AGHS4 0 . . . .

Type III Tests of Fixed Effects

Type III Tests of Fixed Effects
Effect Num DF Den DF F Value Pr > F
Year 4 144 4058.38 <.0001
Location 1 36 90.58 <.0001
Year*Location 4 144 19.76 <.0001
Family 5 36 18.48 <.0001
Year*Family 20 144 4.79 <.0001
Location*Family 5 36 0.13 0.9849
Year*Location*Family 20 144 0.77 0.7482

Solution for Random Effects

Solution for Random Effects
Effect Subject Estimate Std Err Pred DF t Value Pr > |t| Alpha Lower Upper
Intercept GLB AGH25_Hill Farm_1 0.2413 0.2963 46.37 0.81 0.4195 0.05 -0.3549 0.8375
Intercept GLB AGH25_Hill Farm_2 -0.2756 0.2963 46.37 -0.93 0.3570 0.05 -0.8718 0.3206
Intercept GLB AGH25_Hill Farm_3 -0.01982 0.2963 46.37 -0.07 0.9470 0.05 -0.6160 0.5764
Intercept GLB AGH25_Hill Farm_4 0.05412 0.2963 46.37 0.18 0.8558 0.05 -0.5421 0.6503
Intercept GLB AGH25_LA Tech_1 -0.4092 0.2963 46.37 -1.38 0.1738 0.05 -1.0054 0.1870
Intercept GLB AGH25_LA Tech_2 -0.3500 0.2963 46.37 -1.18 0.2434 0.05 -0.9463 0.2462
Intercept GLB AGH25_LA Tech_3 0.3857 0.2963 46.37 1.30 0.1994 0.05 -0.2105 0.9819
Intercept GLB AGH25_LA Tech_4 0.3735 0.2963 46.37 1.26 0.2137 0.05 -0.2227 0.9698
Intercept GLB AGH2_Hill Farm_1 -0.2716 0.2963 46.37 -0.92 0.3640 0.05 -0.8678 0.3246
Intercept GLB AGH2_Hill Farm_2 0.1398 0.2963 46.37 0.47 0.6392 0.05 -0.4564 0.7360
Intercept GLB AGH2_Hill Farm_3 0.3201 0.2963 46.37 1.08 0.2855 0.05 -0.2761 0.9164
Intercept GLB AGH2_Hill Farm_4 -0.1884 0.2963 46.37 -0.64 0.5280 0.05 -0.7846 0.4079
Intercept GLB AGH2_LA Tech_1 -0.05234 0.2963 46.37 -0.18 0.8606 0.05 -0.6485 0.5439
Intercept GLB AGH2_LA Tech_2 -0.6627 0.2963 46.37 -2.24 0.0301 0.05 -1.2590 -0.06653
Intercept GLB AGH2_LA Tech_3 0.8727 0.2963 46.37 2.95 0.0050 0.05 0.2764 1.4689
Intercept GLB AGH2_LA Tech_4 -0.1576 0.2963 46.37 -0.53 0.5973 0.05 -0.7538 0.4386
Intercept GLB AGHS1_Hill Farm_1 0.1936 0.2963 46.37 0.65 0.5166 0.05 -0.4026 0.7898
Intercept GLB AGHS1_Hill Farm_2 -0.1444 0.2963 46.37 -0.49 0.6283 0.05 -0.7406 0.4518
Intercept GLB AGHS1_Hill Farm_3 0.4937 0.2963 46.37 1.67 0.1024 0.05 -0.1025 1.0899
Intercept GLB AGHS1_Hill Farm_4 -0.5429 0.2963 46.37 -1.83 0.0733 0.05 -1.1391 0.05331
Intercept GLB AGHS1_LA Tech_1 -0.4634 0.2963 46.37 -1.56 0.1246 0.05 -1.0596 0.1328
Intercept GLB AGHS1_LA Tech_2 0.4700 0.2963 46.37 1.59 0.1195 0.05 -0.1262 1.0662
Intercept GLB AGHS1_LA Tech_3 0.4483 0.2963 46.37 1.51 0.1370 0.05 -0.1479 1.0445
Intercept GLB AGHS1_LA Tech_4 -0.4549 0.2963 46.37 -1.54 0.1315 0.05 -1.0511 0.1413
Intercept GLB AGHS2_Hill Farm_1 -0.00318 0.2963 46.37 -0.01 0.9915 0.05 -0.5994 0.5930
Intercept GLB AGHS2_Hill Farm_2 -0.8445 0.2963 46.37 -2.85 0.0065 0.05 -1.4407 -0.2483
Intercept GLB AGHS2_Hill Farm_3 0.3695 0.2963 46.37 1.25 0.2185 0.05 -0.2267 0.9657
Intercept GLB AGHS2_Hill Farm_4 0.4782 0.2963 46.37 1.61 0.1133 0.05 -0.1180 1.0744
Intercept GLB AGHS2_LA Tech_1 -0.5213 0.2963 46.37 -1.76 0.0851 0.05 -1.1175 0.07492
Intercept GLB AGHS2_LA Tech_2 0.05283 0.2963 46.37 0.18 0.8592 0.05 -0.5434 0.6490
Intercept GLB AGHS2_LA Tech_3 0.5827 0.2963 46.37 1.97 0.0552 0.05 -0.01349 1.1789
Intercept GLB AGHS2_LA Tech_4 -0.1143 0.2963 46.37 -0.39 0.7015 0.05 -0.7105 0.4820
Intercept GLB AGHS3_Hill Farm_1 0.8690 0.2963 46.37 2.93 0.0052 0.05 0.2727 1.4652
Intercept GLB AGHS3_Hill Farm_2 -1.1000 0.2963 46.37 -3.71 0.0005 0.05 -1.6962 -0.5038
Intercept GLB AGHS3_Hill Farm_3 0.5274 0.2963 46.37 1.78 0.0816 0.05 -0.06881 1.1236
Intercept GLB AGHS3_Hill Farm_4 -0.2963 0.2963 46.37 -1.00 0.3224 0.05 -0.8926 0.2999
Intercept GLB AGHS3_LA Tech_1 -0.2220 0.2963 46.37 -0.75 0.4574 0.05 -0.8182 0.3742
Intercept GLB AGHS3_LA Tech_2 -0.06415 0.2963 46.37 -0.22 0.8295 0.05 -0.6604 0.5321
Intercept GLB AGHS3_LA Tech_3 0.4623 0.2963 46.37 1.56 0.1254 0.05 -0.1339 1.0585
Intercept GLB AGHS3_LA Tech_4 -0.1761 0.2963 46.37 -0.59 0.5550 0.05 -0.7724 0.4201
Intercept GLB AGHS4_Hill Farm_1 0.002375 0.2963 46.37 0.01 0.9936 0.05 -0.5938 0.5986
Intercept GLB AGHS4_Hill Farm_2 0.5383 0.2963 46.37 1.82 0.0757 0.05 -0.05786 1.1346
Intercept GLB AGHS4_Hill Farm_3 -0.2894 0.2963 46.37 -0.98 0.3337 0.05 -0.8856 0.3068
Intercept GLB AGHS4_Hill Farm_4 -0.2513 0.2963 46.37 -0.85 0.4006 0.05 -0.8475 0.3449
Intercept GLB AGHS4_LA Tech_1 -0.2461 0.2963 46.37 -0.83 0.4104 0.05 -0.8423 0.3501
Intercept GLB AGHS4_LA Tech_2 -0.2476 0.2963 46.37 -0.84 0.4076 0.05 -0.8438 0.3486
Intercept GLB AGHS4_LA Tech_3 0.1678 0.2963 46.37 0.57 0.5738 0.05 -0.4284 0.7641
Intercept GLB AGHS4_LA Tech_4 0.3258 0.2963 46.37 1.10 0.2771 0.05 -0.2704 0.9221

Slices

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Fifth 1 68.68 112.48 <.0001

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Location _Location Estimate Standard Error DF t Value Pr > |t| Adj P
Year Fifth Hill Farm LA Tech -1.9848 0.1871 68.68 -10.61 <.0001 <.0001

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Location Estimate  
Year Fifth LA Tech 9.1051 A
Year Fifth      
Year Fifth Hill Farm 7.1203 B

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Fourth 1 68.68 84.50 <.0001

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Location _Location Estimate Standard Error DF t Value Pr > |t| Adj P
Year Fourth Hill Farm LA Tech -1.7203 0.1871 68.68 -9.19 <.0001 <.0001

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Location Estimate  
Year Fourth LA Tech 6.3744 A
Year Fourth      
Year Fourth Hill Farm 4.6541 B

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Seventh 1 68.68 80.17 <.0001

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Location _Location Estimate Standard Error DF t Value Pr > |t| Adj P
Year Seventh Hill Farm LA Tech -1.6756 0.1871 68.68 -8.95 <.0001 <.0001

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Location Estimate  
Year Seventh LA Tech 12.2550 A
Year Seventh      
Year Seventh Hill Farm 10.5794 B

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Sixth 1 68.68 59.98 <.0001

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Location _Location Estimate Standard Error DF t Value Pr > |t| Adj P
Year Sixth Hill Farm LA Tech -1.4494 0.1871 68.68 -7.74 <.0001 <.0001

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Location Estimate  
Year Sixth LA Tech 11.0402 A
Year Sixth      
Year Sixth Hill Farm 9.5908 B

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Third 1 68.68 12.82 0.0006

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Location _Location Estimate Standard Error DF t Value Pr > |t| Adj P
Year Third Hill Farm LA Tech -0.6700 0.1871 68.68 -3.58 0.0006 0.0005

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Location Estimate  
Year Third LA Tech 2.9246 A
Year Third      
Year Third Hill Farm 2.2546 B

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Fifth 5 68.68 14.17 <.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Family _Family Estimate Standard Error DF t Value Pr > |t| Adj P
Year Fifth AGH2 AGH25 -0.4203 0.3241 68.68 -1.30 0.1991 0.7865
Year Fifth AGH2 AGHS1 -0.9569 0.3241 68.68 -2.95 0.0043 0.0420
Year Fifth AGH2 AGHS2 -1.1507 0.3241 68.68 -3.55 0.0007 0.0068
Year Fifth AGH2 AGHS3 -1.2211 0.3241 68.68 -3.77 0.0003 0.0032
Year Fifth AGH2 AGHS4 -2.5250 0.3241 68.68 -7.79 <.0001 <.0001
Year Fifth AGH25 AGHS1 -0.5366 0.3241 68.68 -1.66 0.1024 0.5634
Year Fifth AGH25 AGHS2 -0.7304 0.3241 68.68 -2.25 0.0274 0.2204
Year Fifth AGH25 AGHS3 -0.8007 0.3241 68.68 -2.47 0.0160 0.1400
Year Fifth AGH25 AGHS4 -2.1047 0.3241 68.68 -6.49 <.0001 <.0001
Year Fifth AGHS1 AGHS2 -0.1938 0.3241 68.68 -0.60 0.5519 0.9910
Year Fifth AGHS1 AGHS3 -0.2641 0.3241 68.68 -0.81 0.4179 0.9644
Year Fifth AGHS1 AGHS4 -1.5681 0.3241 68.68 -4.84 <.0001 <.0001
Year Fifth AGHS2 AGHS3 -0.07037 0.3241 68.68 -0.22 0.8288 0.9999
Year Fifth AGHS2 AGHS4 -1.3743 0.3241 68.68 -4.24 <.0001 0.0006
Year Fifth AGHS3 AGHS4 -1.3040 0.3241 68.68 -4.02 0.0001 0.0013

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Family Estimate  
Year Fifth AGHS4 9.5920   A
Year Fifth        
Year Fifth AGHS3 8.2881   B
Year Fifth       B
Year Fifth AGHS2 8.2177   B
Year Fifth       B
Year Fifth AGHS1 8.0239   B
Year Fifth       B
Year Fifth AGH25 7.4873 C B
Year Fifth     C  
Year Fifth AGH2 7.0670 C  

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Fourth 5 68.68 16.51 <.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Family _Family Estimate Standard Error DF t Value Pr > |t| Adj P
Year Fourth AGH2 AGH25 -0.3639 0.3241 68.68 -1.12 0.2654 0.8712
Year Fourth AGH2 AGHS1 -1.6194 0.3241 68.68 -5.00 <.0001 <.0001
Year Fourth AGH2 AGHS2 -1.6211 0.3241 68.68 -5.00 <.0001 <.0001
Year Fourth AGH2 AGHS3 -1.1847 0.3241 68.68 -3.65 0.0005 0.0048
Year Fourth AGH2 AGHS4 -2.5590 0.3241 68.68 -7.89 <.0001 <.0001
Year Fourth AGH25 AGHS1 -1.2555 0.3241 68.68 -3.87 0.0002 0.0022
Year Fourth AGH25 AGHS2 -1.2572 0.3241 68.68 -3.88 0.0002 0.0022
Year Fourth AGH25 AGHS3 -0.8208 0.3241 68.68 -2.53 0.0136 0.1217
Year Fourth AGH25 AGHS4 -2.1950 0.3241 68.68 -6.77 <.0001 <.0001
Year Fourth AGHS1 AGHS2 -0.00171 0.3241 68.68 -0.01 0.9958 1.0000
Year Fourth AGHS1 AGHS3 0.4347 0.3241 68.68 1.34 0.1843 0.7617
Year Fourth AGHS1 AGHS4 -0.9396 0.3241 68.68 -2.90 0.0050 0.0487
Year Fourth AGHS2 AGHS3 0.4364 0.3241 68.68 1.35 0.1826 0.7586
Year Fourth AGHS2 AGHS4 -0.9379 0.3241 68.68 -2.89 0.0051 0.0493
Year Fourth AGHS3 AGHS4 -1.3742 0.3241 68.68 -4.24 <.0001 0.0006

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Family Estimate  
Year Fourth AGHS4 6.8485   A
Year Fourth        
Year Fourth AGHS2 5.9107   B
Year Fourth       B
Year Fourth AGHS1 5.9090   B
Year Fourth       B
Year Fourth AGHS3 5.4743 C B
Year Fourth     C  
Year Fourth AGH25 4.6535 C D
Year Fourth       D
Year Fourth AGH2 4.2896   D

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Seventh 5 68.68 13.26 <.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Family _Family Estimate Standard Error DF t Value Pr > |t| Adj P
Year Seventh AGH2 AGH25 -0.02369 0.3241 68.68 -0.07 0.9420 1.0000
Year Seventh AGH2 AGHS1 -0.4087 0.3241 68.68 -1.26 0.2116 0.8056
Year Seventh AGH2 AGHS2 -0.7000 0.3241 68.68 -2.16 0.0343 0.2634
Year Seventh AGH2 AGHS3 -0.5881 0.3241 68.68 -1.81 0.0740 0.4599
Year Seventh AGH2 AGHS4 -2.2643 0.3241 68.68 -6.99 <.0001 <.0001
Year Seventh AGH25 AGHS1 -0.3850 0.3241 68.68 -1.19 0.2390 0.8420
Year Seventh AGH25 AGHS2 -0.6763 0.3241 68.68 -2.09 0.0407 0.3003
Year Seventh AGH25 AGHS3 -0.5644 0.3241 68.68 -1.74 0.0861 0.5070
Year Seventh AGH25 AGHS4 -2.2406 0.3241 68.68 -6.91 <.0001 <.0001
Year Seventh AGHS1 AGHS2 -0.2912 0.3241 68.68 -0.90 0.3721 0.9463
Year Seventh AGHS1 AGHS3 -0.1794 0.3241 68.68 -0.55 0.5818 0.9937
Year Seventh AGHS1 AGHS4 -1.8555 0.3241 68.68 -5.72 <.0001 <.0001
Year Seventh AGHS2 AGHS3 0.1119 0.3241 68.68 0.35 0.7310 0.9993
Year Seventh AGHS2 AGHS4 -1.5643 0.3241 68.68 -4.83 <.0001 <.0001
Year Seventh AGHS3 AGHS4 -1.6762 0.3241 68.68 -5.17 <.0001 <.0001

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Family Estimate  
Year Seventh AGHS4 13.0174 A
Year Seventh      
Year Seventh AGHS2 11.4531 B
Year Seventh     B
Year Seventh AGHS3 11.3412 B
Year Seventh     B
Year Seventh AGHS1 11.1618 B
Year Seventh     B
Year Seventh AGH25 10.7768 B
Year Seventh     B
Year Seventh AGH2 10.7531 B

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Sixth 5 68.68 16.04 <.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Family _Family Estimate Standard Error DF t Value Pr > |t| Adj P
Year Sixth AGH2 AGH25 -0.2837 0.3241 68.68 -0.88 0.3845 0.9519
Year Sixth AGH2 AGHS1 -0.9029 0.3241 68.68 -2.79 0.0069 0.0656
Year Sixth AGH2 AGHS2 -0.9487 0.3241 68.68 -2.93 0.0046 0.0451
Year Sixth AGH2 AGHS3 -1.3610 0.3241 68.68 -4.20 <.0001 0.0007
Year Sixth AGH2 AGHS4 -2.5997 0.3241 68.68 -8.02 <.0001 <.0001
Year Sixth AGH25 AGHS1 -0.6192 0.3241 68.68 -1.91 0.0603 0.4002
Year Sixth AGH25 AGHS2 -0.6650 0.3241 68.68 -2.05 0.0440 0.3187
Year Sixth AGH25 AGHS3 -1.0774 0.3241 68.68 -3.32 0.0014 0.0141
Year Sixth AGH25 AGHS4 -2.3160 0.3241 68.68 -7.15 <.0001 <.0001
Year Sixth AGHS1 AGHS2 -0.04586 0.3241 68.68 -0.14 0.8879 1.0000
Year Sixth AGHS1 AGHS3 -0.4582 0.3241 68.68 -1.41 0.1620 0.7188
Year Sixth AGHS1 AGHS4 -1.6969 0.3241 68.68 -5.23 <.0001 <.0001
Year Sixth AGHS2 AGHS3 -0.4123 0.3241 68.68 -1.27 0.2077 0.7998
Year Sixth AGHS2 AGHS4 -1.6510 0.3241 68.68 -5.09 <.0001 <.0001
Year Sixth AGHS3 AGHS4 -1.2387 0.3241 68.68 -3.82 0.0003 0.0027

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Family Estimate  
Year Sixth AGHS4 11.8992   A  
Year Sixth          
Year Sixth AGHS3 10.6605   B  
Year Sixth       B  
Year Sixth AGHS2 10.2482 C B  
Year Sixth     C B  
Year Sixth AGHS1 10.2023 C B D
Year Sixth     C   D
Year Sixth AGH25 9.5831 C   D
Year Sixth         D
Year Sixth AGH2 9.2995     D

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Third 5 68.68 12.50 <.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Family _Family Estimate Standard Error DF t Value Pr > |t| Adj P
Year Third AGH2 AGH25 -0.09482 0.3241 68.68 -0.29 0.7708 0.9997
Year Third AGH2 AGHS1 -1.5896 0.3241 68.68 -4.90 <.0001 <.0001
Year Third AGH2 AGHS2 -1.6557 0.3241 68.68 -5.11 <.0001 <.0001
Year Third AGH2 AGHS3 -0.6881 0.3241 68.68 -2.12 0.0374 0.2815
Year Third AGH2 AGHS4 -1.7722 0.3241 68.68 -5.47 <.0001 <.0001
Year Third AGH25 AGHS1 -1.4948 0.3241 68.68 -4.61 <.0001 0.0001
Year Third AGH25 AGHS2 -1.5609 0.3241 68.68 -4.82 <.0001 <.0001
Year Third AGH25 AGHS3 -0.5933 0.3241 68.68 -1.83 0.0715 0.4497
Year Third AGH25 AGHS4 -1.6774 0.3241 68.68 -5.17 <.0001 <.0001
Year Third AGHS1 AGHS2 -0.06610 0.3241 68.68 -0.20 0.8390 1.0000
Year Third AGHS1 AGHS3 0.9015 0.3241 68.68 2.78 0.0070 0.0663
Year Third AGHS1 AGHS4 -0.1826 0.3241 68.68 -0.56 0.5751 0.9932
Year Third AGHS2 AGHS3 0.9676 0.3241 68.68 2.99 0.0039 0.0384
Year Third AGHS2 AGHS4 -0.1165 0.3241 68.68 -0.36 0.7204 0.9992
Year Third AGHS3 AGHS4 -1.0841 0.3241 68.68 -3.34 0.0013 0.0132

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Family Estimate  
Year Third AGHS4 3.3951   A
Year Third       A
Year Third AGHS2 3.2786   A
Year Third       A
Year Third AGHS1 3.2125 B A
Year Third     B  
Year Third AGHS3 2.3110 B C
Year Third       C
Year Third AGH25 1.7177   C
Year Third       C
Year Third AGH2 1.6229   C

LS-Means Results

Family LS-Means

Family Least Squares Means
Family Estimate Standard
Error
DF t Value Pr > |t|
AGH2 6.6064 0.1930 36 34.23 <.0001
AGH25 6.8437 0.1930 36 35.45 <.0001
AGHS1 7.7019 0.1930 36 39.90 <.0001
AGHS2 7.8216 0.1930 36 40.52 <.0001
AGHS3 7.6150 0.1930 36 39.45 <.0001
AGHS4 8.9504 0.1930 36 46.37 <.0001

Family Diffs

Differences of Family Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Family _Family Estimate Standard Error DF t Value Pr > |t| Adj P
AGH2 AGH25 -0.2373 0.2730 36 -0.87 0.3905 0.9514
AGH2 AGHS1 -1.0955 0.2730 36 -4.01 0.0003 0.0037
AGH2 AGHS2 -1.2152 0.2730 36 -4.45 <.0001 0.0010
AGH2 AGHS3 -1.0086 0.2730 36 -3.69 0.0007 0.0088
AGH2 AGHS4 -2.3440 0.2730 36 -8.59 <.0001 <.0001
AGH25 AGHS1 -0.8582 0.2730 36 -3.14 0.0033 0.0362
AGH25 AGHS2 -0.9780 0.2730 36 -3.58 0.0010 0.0118
AGH25 AGHS3 -0.7713 0.2730 36 -2.83 0.0077 0.0760
AGH25 AGHS4 -2.1067 0.2730 36 -7.72 <.0001 <.0001
AGHS1 AGHS2 -0.1197 0.2730 36 -0.44 0.6635 0.9978
AGHS1 AGHS3 0.08690 0.2730 36 0.32 0.7521 0.9995
AGHS1 AGHS4 -1.2485 0.2730 36 -4.57 <.0001 0.0007
AGHS2 AGHS3 0.2066 0.2730 36 0.76 0.4540 0.9729
AGHS2 AGHS4 -1.1288 0.2730 36 -4.14 0.0002 0.0026
AGHS3 AGHS4 -1.3354 0.2730 36 -4.89 <.0001 0.0003

Family Lines

Tukey Grouping for Family Least Squares Means (Alpha=0.05)
LS-means with the same letter are not significantly different.
Family Estimate  
AGHS4 8.9504   A
       
AGHS2 7.8216   B
      B
AGHS1 7.7019   B
      B
AGHS3 7.6150 C B
    C  
AGH25 6.8437 C D
      D
AGH2 6.6064   D

SAS Output Height

htmltools::includeHTML("https://raw.githubusercontent.com/azd169/SAS_outputs/main/hane-results-height.html")
## Warning: `includeHTML()` was provided a `path` that appears to be a complete HTML document.
## ✖ Path: https://raw.githubusercontent.com/azd169/SAS_outputs/main/hane-results-height.html
## ℹ Use `tags$iframe()` to include an HTML document. You can either ensure `path` is accessible in your app or document (see e.g. `shiny::addResourcePath()`) and pass the relative path to the `src` argument. Or you can read the contents of `path` and pass the contents to `srcdoc`.
Results: hane.sas

Results: hane.sas

The GLIMMIX Procedure

The Glimmix Procedure

Model Information

Model Information
Data Set WORK.HEIGHT
Response Variable Height
Response Distribution Gaussian
Link Function Identity
Variance Function Default
Variance Matrix Blocked By GLB
Estimation Technique Restricted Maximum Likelihood
Degrees of Freedom Method Kenward-Roger
Fixed Effects SE Adjustment Kenward-Roger

Class Level Information

Class Level Information
Class Levels Values
Year 8 Fifth First Fourth Initial Second Seventh Sixth Third
Location 2 Hill Farm LA Tech
Block 4 1 2 3 4
Family 6 AGH2 AGH25 AGHS1 AGHS2 AGHS3 AGHS4
GLB 48 AGH25_Hill Farm_1 AGH25_Hill Farm_2 AGH25_Hill Farm_3 AGH25_Hill Farm_4 AGH25_LA Tech_1 AGH25_LA Tech_2 AGH25_LA Tech_3 AGH25_LA Tech_4 AGH2_Hill Farm_1 AGH2_Hill Farm_2 AGH2_Hill Farm_3 AGH2_Hill Farm_4 AGH2_LA Tech_1 AGH2_LA Tech_2 AGH2_LA Tech_3 AGH2_LA Tech_4 AGHS1_Hill Farm_1 AGHS1_Hill Farm_2 AGHS1_Hill Farm_3 AGHS1_Hill Farm_4 AGHS1_LA Tech_1 AGHS1_LA Tech_2 AGHS1_LA Tech_3 AGHS1_LA Tech_4 AGHS2_Hill Farm_1 AGHS2_Hill Farm_2 AGHS2_Hill Farm_3 AGHS2_Hill Farm_4 AGHS2_LA Tech_1 AGHS2_LA Tech_2 AGHS2_LA Tech_3 AGHS2_LA Tech_4 AGHS3_Hill Farm_1 AGHS3_Hill Farm_2 AGHS3_Hill Farm_3 AGHS3_Hill Farm_4 AGHS3_LA Tech_1 AGHS3_LA Tech_2 AGHS3_LA Tech_3 AGHS3_LA Tech_4 AGHS4_Hill Farm_1 AGHS4_Hill Farm_2 AGHS4_Hill Farm_3 AGHS4_Hill Farm_4 AGHS4_LA Tech_1 AGHS4_LA Tech_2 AGHS4_LA Tech_3 AGHS4_LA Tech_4

Number of Observations

Number of Observations Read 384
Number of Observations Used 384

Dimensions

Dimensions
G-side Cov. Parameters 2
R-side Cov. Parameters 1
Columns in X 189
Columns in Z per Subject 1
Subjects (Blocks in V) 48
Max Obs per Subject 8

Optimization Information

Optimization Information
Optimization Technique Dual Quasi-Newton
Parameters in Optimization 2
Lower Boundaries 2
Upper Boundaries 1
Fixed Effects Profiled
Residual Variance Profiled
Starting From Data

Iteration History

Iteration History
Iteration Restarts Evaluations Objective
Function
Change Max
Gradient
0 0 4 419.85607123 . 3.56E-12

Convergence Status

Convergence criterion (ABSGCONV=0.00001) satisfied.

Hessian Matrix

Hessian Matrix
Cov Parm Subject Row Col1 Col2
Variance GLB 1 80.0040  
AR(1) GLB 2    

Fit Statistics

Fit Statistics
-2 Res Log Likelihood 419.86
AIC (smaller is better) 425.86
AICC (smaller is better) 425.94
BIC (smaller is better) 431.47
CAIC (smaller is better) 434.47
HQIC (smaller is better) 427.98
Generalized Chi-Square 37.31
Gener. Chi-Square / DF 0.13

Covariance Parameter Estimates

Covariance Parameter Estimates
Cov Parm Subject Estimate Standard
Error
Variance GLB 0.06509 0.01921
AR(1) GLB 0 .
Residual   0.1295 0.01154

Solutions for Fixed Effects

Solutions for Fixed Effects
Effect Year Location Family Estimate Standard
Error
DF t Value Pr > |t|
Intercept       3.6033 0.2206 161.5 16.34 <.0001
Year Fifth     3.7856 0.2545 252 14.88 <.0001
Year First     -2.6121 0.2545 252 -10.26 <.0001
Year Fourth     2.2226 0.2545 252 8.73 <.0001
Year Initial     -3.0528 0.2545 252 -12.00 <.0001
Year Second     -1.2216 0.2545 252 -4.80 <.0001
Year Seventh     5.8582 0.2545 252 23.02 <.0001
Year Sixth     5.9387 0.2545 252 23.34 <.0001
Year Third     0 . . . .
Location   Hill Farm   -0.7677 0.3119 161.5 -2.46 0.0149
Location   LA Tech   0 . . . .
Year*Location Fifth Hill Farm   -0.3317 0.3599 252 -0.92 0.3575
Year*Location Fifth LA Tech   0 . . . .
Year*Location First Hill Farm   0.7031 0.3599 252 1.95 0.0518
Year*Location First LA Tech   0 . . . .
Year*Location Fourth Hill Farm   -0.5411 0.3599 252 -1.50 0.1340
Year*Location Fourth LA Tech   0 . . . .
Year*Location Initial Hill Farm   0.7707 0.3599 252 2.14 0.0332
Year*Location Initial LA Tech   0 . . . .
Year*Location Second Hill Farm   0.4508 0.3599 252 1.25 0.2115
Year*Location Second LA Tech   0 . . . .
Year*Location Seventh Hill Farm   1.0020 0.3599 252 2.78 0.0058
Year*Location Seventh LA Tech   0 . . . .
Year*Location Sixth Hill Farm   -2.5706 0.3599 252 -7.14 <.0001
Year*Location Sixth LA Tech   0 . . . .
Year*Location Third Hill Farm   0 . . . .
Year*Location Third LA Tech   0 . . . .
Family     AGH2 -0.9547 0.3119 161.5 -3.06 0.0026
Family     AGH25 -0.7602 0.3119 161.5 -2.44 0.0159
Family     AGHS1 0.2949 0.3119 161.5 0.95 0.3459
Family     AGHS2 0.4191 0.3119 161.5 1.34 0.1810
Family     AGHS3 -0.3812 0.3119 161.5 -1.22 0.2235
Family     AGHS4 0 . . . .
Year*Family Fifth   AGH2 -0.4696 0.3599 252 -1.30 0.1931
Year*Family Fifth   AGH25 -0.04786 0.3599 252 -0.13 0.8943
Year*Family Fifth   AGHS1 -0.1533 0.3599 252 -0.43 0.6706
Year*Family Fifth   AGHS2 -0.05521 0.3599 252 -0.15 0.8782
Year*Family Fifth   AGHS3 -0.09546 0.3599 252 -0.27 0.7910
Year*Family Fifth   AGHS4 0 . . . .
Year*Family First   AGH2 0.6650 0.3599 252 1.85 0.0658
Year*Family First   AGH25 0.4395 0.3599 252 1.22 0.2232
Year*Family First   AGHS1 -0.3370 0.3599 252 -0.94 0.3500
Year*Family First   AGHS2 -0.4107 0.3599 252 -1.14 0.2549
Year*Family First   AGHS3 0.2685 0.3599 252 0.75 0.4563
Year*Family First   AGHS4 0 . . . .
Year*Family Fourth   AGH2 -0.5640 0.3599 252 -1.57 0.1184
Year*Family Fourth   AGH25 -0.3129 0.3599 252 -0.87 0.3855
Year*Family Fourth   AGHS1 -0.06452 0.3599 252 -0.18 0.8579
Year*Family Fourth   AGHS2 0.07452 0.3599 252 0.21 0.8361
Year*Family Fourth   AGHS3 0.06692 0.3599 252 0.19 0.8526
Year*Family Fourth   AGHS4 0 . . . .
Year*Family Initial   AGH2 0.8553 0.3599 252 2.38 0.0182
Year*Family Initial   AGH25 0.6140 0.3599 252 1.71 0.0892
Year*Family Initial   AGHS1 -0.3388 0.3599 252 -0.94 0.3474
Year*Family Initial   AGHS2 -0.5205 0.3599 252 -1.45 0.1494
Year*Family Initial   AGHS3 0.3439 0.3599 252 0.96 0.3402
Year*Family Initial   AGHS4 0 . . . .
Year*Family Second   AGH2 0.1622 0.3599 252 0.45 0.6525
Year*Family Second   AGH25 0.1249 0.3599 252 0.35 0.7288
Year*Family Second   AGHS1 -0.1838 0.3599 252 -0.51 0.6100
Year*Family Second   AGHS2 -0.05966 0.3599 252 -0.17 0.8685
Year*Family Second   AGHS3 0.1786 0.3599 252 0.50 0.6201
Year*Family Second   AGHS4 0 . . . .
Year*Family Seventh   AGH2 -0.1003 0.3599 252 -0.28 0.7808
Year*Family Seventh   AGH25 -0.1628 0.3599 252 -0.45 0.6514
Year*Family Seventh   AGHS1 -0.2814 0.3599 252 -0.78 0.4351
Year*Family Seventh   AGHS2 -0.3133 0.3599 252 -0.87 0.3849
Year*Family Seventh   AGHS3 0.3029 0.3599 252 0.84 0.4009
Year*Family Seventh   AGHS4 0 . . . .
Year*Family Sixth   AGH2 -1.4486 0.3599 252 -4.02 <.0001
Year*Family Sixth   AGH25 -0.9490 0.3599 252 -2.64 0.0089
Year*Family Sixth   AGHS1 -0.06103 0.3599 252 -0.17 0.8655
Year*Family Sixth   AGHS2 -0.1143 0.3599 252 -0.32 0.7511
Year*Family Sixth   AGHS3 -0.05486 0.3599 252 -0.15 0.8790
Year*Family Sixth   AGHS4 0 . . . .
Year*Family Third   AGH2 0 . . . .
Year*Family Third   AGH25 0 . . . .
Year*Family Third   AGHS1 0 . . . .
Year*Family Third   AGHS2 0 . . . .
Year*Family Third   AGHS3 0 . . . .
Year*Family Third   AGHS4 0 . . . .
Location*Family   Hill Farm AGH2 0.4218 0.4412 161.5 0.96 0.3404
Location*Family   Hill Farm AGH25 0.3009 0.4412 161.5 0.68 0.4962
Location*Family   Hill Farm AGHS1 -0.02870 0.4412 161.5 -0.07 0.9482
Location*Family   Hill Farm AGHS2 -0.3231 0.4412 161.5 -0.73 0.4650
Location*Family   Hill Farm AGHS3 -0.05930 0.4412 161.5 -0.13 0.8932
Location*Family   Hill Farm AGHS4 0 . . . .
Location*Family   LA Tech AGH2 0 . . . .
Location*Family   LA Tech AGH25 0 . . . .
Location*Family   LA Tech AGHS1 0 . . . .
Location*Family   LA Tech AGHS2 0 . . . .
Location*Family   LA Tech AGHS3 0 . . . .
Location*Family   LA Tech AGHS4 0 . . . .
Year*Location*Family Fifth Hill Farm AGH2 -0.4498 0.5090 252 -0.88 0.3776
Year*Location*Family Fifth Hill Farm AGH25 -0.5438 0.5090 252 -1.07 0.2863
Year*Location*Family Fifth Hill Farm AGHS1 -0.3138 0.5090 252 -0.62 0.5381
Year*Location*Family Fifth Hill Farm AGHS2 0.1410 0.5090 252 0.28 0.7820
Year*Location*Family Fifth Hill Farm AGHS3 -0.3825 0.5090 252 -0.75 0.4531
Year*Location*Family Fifth Hill Farm AGHS4 0 . . . .
Year*Location*Family Fifth LA Tech AGH2 0 . . . .
Year*Location*Family Fifth LA Tech AGH25 0 . . . .
Year*Location*Family Fifth LA Tech AGHS1 0 . . . .
Year*Location*Family Fifth LA Tech AGHS2 0 . . . .
Year*Location*Family Fifth LA Tech AGHS3 0 . . . .
Year*Location*Family Fifth LA Tech AGHS4 0 . . . .
Year*Location*Family First Hill Farm AGH2 -0.3662 0.5090 252 -0.72 0.4725
Year*Location*Family First Hill Farm AGH25 -0.1775 0.5090 252 -0.35 0.7276
Year*Location*Family First Hill Farm AGHS1 0.1585 0.5090 252 0.31 0.7557
Year*Location*Family First Hill Farm AGHS2 0.3192 0.5090 252 0.63 0.5312
Year*Location*Family First Hill Farm AGHS3 0.1219 0.5090 252 0.24 0.8109
Year*Location*Family First Hill Farm AGHS4 0 . . . .
Year*Location*Family First LA Tech AGH2 0 . . . .
Year*Location*Family First LA Tech AGH25 0 . . . .
Year*Location*Family First LA Tech AGHS1 0 . . . .
Year*Location*Family First LA Tech AGHS2 0 . . . .
Year*Location*Family First LA Tech AGHS3 0 . . . .
Year*Location*Family First LA Tech AGHS4 0 . . . .
Year*Location*Family Fourth Hill Farm AGH2 -0.1764 0.5090 252 -0.35 0.7292
Year*Location*Family Fourth Hill Farm AGH25 -0.2610 0.5090 252 -0.51 0.6086
Year*Location*Family Fourth Hill Farm AGHS1 -0.08584 0.5090 252 -0.17 0.8662
Year*Location*Family Fourth Hill Farm AGHS2 0.07420 0.5090 252 0.15 0.8842
Year*Location*Family Fourth Hill Farm AGHS3 -0.2134 0.5090 252 -0.42 0.6754
Year*Location*Family Fourth Hill Farm AGHS4 0 . . . .
Year*Location*Family Fourth LA Tech AGH2 0 . . . .
Year*Location*Family Fourth LA Tech AGH25 0 . . . .
Year*Location*Family Fourth LA Tech AGHS1 0 . . . .
Year*Location*Family Fourth LA Tech AGHS2 0 . . . .
Year*Location*Family Fourth LA Tech AGHS3 0 . . . .
Year*Location*Family Fourth LA Tech AGHS4 0 . . . .
Year*Location*Family Initial Hill Farm AGH2 -0.4515 0.5090 252 -0.89 0.3759
Year*Location*Family Initial Hill Farm AGH25 -0.2910 0.5090 252 -0.57 0.5681
Year*Location*Family Initial Hill Farm AGHS1 0.08518 0.5090 252 0.17 0.8672
Year*Location*Family Initial Hill Farm AGHS2 0.3375 0.5090 252 0.66 0.5079
Year*Location*Family Initial Hill Farm AGHS3 0.08876 0.5090 252 0.17 0.8617
Year*Location*Family Initial Hill Farm AGHS4 0 . . . .
Year*Location*Family Initial LA Tech AGH2 0 . . . .
Year*Location*Family Initial LA Tech AGH25 0 . . . .
Year*Location*Family Initial LA Tech AGHS1 0 . . . .
Year*Location*Family Initial LA Tech AGHS2 0 . . . .
Year*Location*Family Initial LA Tech AGHS3 0 . . . .
Year*Location*Family Initial LA Tech AGHS4 0 . . . .
Year*Location*Family Second Hill Farm AGH2 -0.1221 0.5090 252 -0.24 0.8106
Year*Location*Family Second Hill Farm AGH25 -0.1298 0.5090 252 -0.26 0.7989
Year*Location*Family Second Hill Farm AGHS1 0.06376 0.5090 252 0.13 0.9004
Year*Location*Family Second Hill Farm AGHS2 0.1648 0.5090 252 0.32 0.7463
Year*Location*Family Second Hill Farm AGHS3 0.1828 0.5090 252 0.36 0.7197
Year*Location*Family Second Hill Farm AGHS4 0 . . . .
Year*Location*Family Second LA Tech AGH2 0 . . . .
Year*Location*Family Second LA Tech AGH25 0 . . . .
Year*Location*Family Second LA Tech AGHS1 0 . . . .
Year*Location*Family Second LA Tech AGHS2 0 . . . .
Year*Location*Family Second LA Tech AGHS3 0 . . . .
Year*Location*Family Second LA Tech AGHS4 0 . . . .
Year*Location*Family Seventh Hill Farm AGH2 -0.6213 0.5090 252 -1.22 0.2234
Year*Location*Family Seventh Hill Farm AGH25 -0.9667 0.5090 252 -1.90 0.0587
Year*Location*Family Seventh Hill Farm AGHS1 0.1514 0.5090 252 0.30 0.7664
Year*Location*Family Seventh Hill Farm AGHS2 0.2505 0.5090 252 0.49 0.6230
Year*Location*Family Seventh Hill Farm AGHS3 -1.0173 0.5090 252 -2.00 0.0467
Year*Location*Family Seventh Hill Farm AGHS4 0 . . . .
Year*Location*Family Seventh LA Tech AGH2 0 . . . .
Year*Location*Family Seventh LA Tech AGH25 0 . . . .
Year*Location*Family Seventh LA Tech AGHS1 0 . . . .
Year*Location*Family Seventh LA Tech AGHS2 0 . . . .
Year*Location*Family Seventh LA Tech AGHS3 0 . . . .
Year*Location*Family Seventh LA Tech AGHS4 0 . . . .
Year*Location*Family Sixth Hill Farm AGH2 1.1240 0.5090 252 2.21 0.0281
Year*Location*Family Sixth Hill Farm AGH25 0.5276 0.5090 252 1.04 0.3010
Year*Location*Family Sixth Hill Farm AGHS1 0.06816 0.5090 252 0.13 0.8936
Year*Location*Family Sixth Hill Farm AGHS2 0.3570 0.5090 252 0.70 0.4838
Year*Location*Family Sixth Hill Farm AGHS3 0.1141 0.5090 252 0.22 0.8228
Year*Location*Family Sixth Hill Farm AGHS4 0 . . . .
Year*Location*Family Sixth LA Tech AGH2 0 . . . .
Year*Location*Family Sixth LA Tech AGH25 0 . . . .
Year*Location*Family Sixth LA Tech AGHS1 0 . . . .
Year*Location*Family Sixth LA Tech AGHS2 0 . . . .
Year*Location*Family Sixth LA Tech AGHS3 0 . . . .
Year*Location*Family Sixth LA Tech AGHS4 0 . . . .
Year*Location*Family Third Hill Farm AGH2 0 . . . .
Year*Location*Family Third Hill Farm AGH25 0 . . . .
Year*Location*Family Third Hill Farm AGHS1 0 . . . .
Year*Location*Family Third Hill Farm AGHS2 0 . . . .
Year*Location*Family Third Hill Farm AGHS3 0 . . . .
Year*Location*Family Third Hill Farm AGHS4 0 . . . .
Year*Location*Family Third LA Tech AGH2 0 . . . .
Year*Location*Family Third LA Tech AGH25 0 . . . .
Year*Location*Family Third LA Tech AGHS1 0 . . . .
Year*Location*Family Third LA Tech AGHS2 0 . . . .
Year*Location*Family Third LA Tech AGHS3 0 . . . .
Year*Location*Family Third LA Tech AGHS4 0 . . . .

Type III Tests of Fixed Effects

Type III Tests of Fixed Effects
Effect Num DF Den DF F Value Pr > F
Year 7 252 3715.43 <.0001
Location 1 36 100.98 <.0001
Year*Location 7 252 95.45 <.0001
Family 5 36 21.61 <.0001
Year*Family 35 252 5.27 <.0001
Location*Family 5 36 0.70 0.6304
Year*Location*Family 35 252 1.24 0.1731

Solution for Random Effects

Solution for Random Effects
Effect Subject Estimate Std Err Pred DF t Value Pr > |t| Alpha Lower Upper
Intercept GLB AGH25_Hill Farm_1 -0.1321 0.1622 46.92 -0.81 0.4196 0.05 -0.4584 0.1942
Intercept GLB AGH25_Hill Farm_2 0.06641 0.1622 46.92 0.41 0.6841 0.05 -0.2599 0.3927
Intercept GLB AGH25_Hill Farm_3 -0.09112 0.1622 46.92 -0.56 0.5769 0.05 -0.4174 0.2352
Intercept GLB AGH25_Hill Farm_4 0.1568 0.1622 46.92 0.97 0.3387 0.05 -0.1695 0.4831
Intercept GLB AGH25_LA Tech_1 -0.1835 0.1622 46.92 -1.13 0.2637 0.05 -0.5098 0.1428
Intercept GLB AGH25_LA Tech_2 -0.1597 0.1622 46.92 -0.98 0.3299 0.05 -0.4860 0.1666
Intercept GLB AGH25_LA Tech_3 0.07326 0.1622 46.92 0.45 0.6536 0.05 -0.2530 0.3995
Intercept GLB AGH25_LA Tech_4 0.2699 0.1622 46.92 1.66 0.1027 0.05 -0.05639 0.5962
Intercept GLB AGH2_Hill Farm_1 -0.07375 0.1622 46.92 -0.45 0.6514 0.05 -0.4000 0.2525
Intercept GLB AGH2_Hill Farm_2 -0.08202 0.1622 46.92 -0.51 0.6154 0.05 -0.4083 0.2443
Intercept GLB AGH2_Hill Farm_3 0.1422 0.1622 46.92 0.88 0.3849 0.05 -0.1841 0.4685
Intercept GLB AGH2_Hill Farm_4 0.01353 0.1622 46.92 0.08 0.9339 0.05 -0.3128 0.3398
Intercept GLB AGH2_LA Tech_1 -0.1325 0.1622 46.92 -0.82 0.4180 0.05 -0.4588 0.1938
Intercept GLB AGH2_LA Tech_2 -0.1581 0.1622 46.92 -0.97 0.3346 0.05 -0.4844 0.1682
Intercept GLB AGH2_LA Tech_3 0.4081 0.1622 46.92 2.52 0.0153 0.05 0.08183 0.7344
Intercept GLB AGH2_LA Tech_4 -0.1175 0.1622 46.92 -0.72 0.4724 0.05 -0.4438 0.2088
Intercept GLB AGHS1_Hill Farm_1 0.03870 0.1622 46.92 0.24 0.8124 0.05 -0.2876 0.3650
Intercept GLB AGHS1_Hill Farm_2 0.03476 0.1622 46.92 0.21 0.8312 0.05 -0.2915 0.3611
Intercept GLB AGHS1_Hill Farm_3 0.2445 0.1622 46.92 1.51 0.1384 0.05 -0.08179 0.5708
Intercept GLB AGHS1_Hill Farm_4 -0.3180 0.1622 46.92 -1.96 0.0559 0.05 -0.6443 0.008324
Intercept GLB AGHS1_LA Tech_1 -0.2625 0.1622 46.92 -1.62 0.1123 0.05 -0.5888 0.06382
Intercept GLB AGHS1_LA Tech_2 0.03883 0.1622 46.92 0.24 0.8118 0.05 -0.2875 0.3651
Intercept GLB AGHS1_LA Tech_3 0.3097 0.1622 46.92 1.91 0.0623 0.05 -0.01662 0.6360
Intercept GLB AGHS1_LA Tech_4 -0.08603 0.1622 46.92 -0.53 0.5983 0.05 -0.4123 0.2403
Intercept GLB AGHS2_Hill Farm_1 -0.2649 0.1622 46.92 -1.63 0.1091 0.05 -0.5912 0.06140
Intercept GLB AGHS2_Hill Farm_2 -0.2065 0.1622 46.92 -1.27 0.2093 0.05 -0.5328 0.1198
Intercept GLB AGHS2_Hill Farm_3 0.2849 0.1622 46.92 1.76 0.0856 0.05 -0.04144 0.6111
Intercept GLB AGHS2_Hill Farm_4 0.1865 0.1622 46.92 1.15 0.2560 0.05 -0.1398 0.5128
Intercept GLB AGHS2_LA Tech_1 -0.2952 0.1622 46.92 -1.82 0.0751 0.05 -0.6215 0.03107
Intercept GLB AGHS2_LA Tech_2 -0.01269 0.1622 46.92 -0.08 0.9380 0.05 -0.3390 0.3136
Intercept GLB AGHS2_LA Tech_3 0.3123 0.1622 46.92 1.93 0.0602 0.05 -0.01401 0.6386
Intercept GLB AGHS2_LA Tech_4 -0.00437 0.1622 46.92 -0.03 0.9786 0.05 -0.3307 0.3219
Intercept GLB AGHS3_Hill Farm_1 -0.00285 0.1622 46.92 -0.02 0.9861 0.05 -0.3291 0.3234
Intercept GLB AGHS3_Hill Farm_2 -0.4217 0.1622 46.92 -2.60 0.0124 0.05 -0.7480 -0.09539
Intercept GLB AGHS3_Hill Farm_3 0.3871 0.1622 46.92 2.39 0.0211 0.05 0.06081 0.7134
Intercept GLB AGHS3_Hill Farm_4 0.03743 0.1622 46.92 0.23 0.8185 0.05 -0.2889 0.3637
Intercept GLB AGHS3_LA Tech_1 -0.2095 0.1622 46.92 -1.29 0.2029 0.05 -0.5357 0.1168
Intercept GLB AGHS3_LA Tech_2 -0.09466 0.1622 46.92 -0.58 0.5623 0.05 -0.4209 0.2316
Intercept GLB AGHS3_LA Tech_3 0.1860 0.1622 46.92 1.15 0.2573 0.05 -0.1403 0.5123
Intercept GLB AGHS3_LA Tech_4 0.1181 0.1622 46.92 0.73 0.4701 0.05 -0.2082 0.4444
Intercept GLB AGHS4_Hill Farm_1 -0.2830 0.1622 46.92 -1.75 0.0875 0.05 -0.6093 0.04324
Intercept GLB AGHS4_Hill Farm_2 0.1853 0.1622 46.92 1.14 0.2590 0.05 -0.1410 0.5116
Intercept GLB AGHS4_Hill Farm_3 -0.00433 0.1622 46.92 -0.03 0.9788 0.05 -0.3306 0.3220
Intercept GLB AGHS4_Hill Farm_4 0.1021 0.1622 46.92 0.63 0.5322 0.05 -0.2242 0.4284
Intercept GLB AGHS4_LA Tech_1 -0.2516 0.1622 46.92 -1.55 0.1275 0.05 -0.5779 0.07465
Intercept GLB AGHS4_LA Tech_2 -0.07550 0.1622 46.92 -0.47 0.6437 0.05 -0.4018 0.2508
Intercept GLB AGHS4_LA Tech_3 0.1680 0.1622 46.92 1.04 0.3055 0.05 -0.1583 0.4943
Intercept GLB AGHS4_LA Tech_4 0.1591 0.1622 46.92 0.98 0.3316 0.05 -0.1672 0.4854

Slices

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Fifth 1 161.5 105.11 <.0001

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Location _Location Estimate Standard Error DF t Value Pr > |t| Adj P
Year Fifth Hill Farm LA Tech -1.3057 0.1274 161.5 -10.25 <.0001 <.0001

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Location Estimate  
Year Fifth LA Tech 7.0217 A
Year Fifth      
Year Fifth Hill Farm 5.7160 B

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year First 1 161.5 0.00 0.9793

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Location _Location Estimate Standard Error DF t Value Pr > |t| Adj P
Year First Hill Farm LA Tech -0.00330 0.1274 161.5 -0.03 0.9793 0.9793

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Location Estimate  
Year First LA Tech 0.8650 A
Year First     A
Year First Hill Farm 0.8617 A

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Fourth 1 161.5 115.26 <.0001

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Location _Location Estimate Standard Error DF t Value Pr > |t| Adj P
Year Fourth Hill Farm LA Tech -1.3672 0.1274 161.5 -10.74 <.0001 <.0001

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Location Estimate  
Year Fourth LA Tech 5.4623 A
Year Fourth      
Year Fourth Hill Farm 4.0950 B

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Initial 1 161.5 0.02 0.8976

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Location _Location Estimate Standard Error DF t Value Pr > |t| Adj P
Year Initial Hill Farm LA Tech 0.01642 0.1274 161.5 0.13 0.8976 0.8975

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Location Estimate  
Year Initial Hill Farm 0.4955 A
Year Initial     A
Year Initial LA Tech 0.4791 A

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Second 1 161.5 3.50 0.0631

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Location _Location Estimate Standard Error DF t Value Pr > |t| Adj P
Year Second Hill Farm LA Tech -0.2384 0.1274 161.5 -1.87 0.0631 0.0624

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Location Estimate  
Year Second LA Tech 2.1884 A
Year Second     A
Year Second Hill Farm 1.9500 A

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Seventh 1 161.5 0.40 0.5259

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Location _Location Estimate Standard Error DF t Value Pr > |t| Adj P
Year Seventh Hill Farm LA Tech -0.08095 0.1274 161.5 -0.64 0.5259 0.5256

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Location Estimate  
Year Seventh LA Tech 9.1387 A
Year Seventh     A
Year Seventh Hill Farm 9.0577 A

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Sixth 1 161.5 526.16 <.0001

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Location _Location Estimate Standard Error DF t Value Pr > |t| Adj P
Year Sixth Hill Farm LA Tech -2.9212 0.1274 161.5 -22.94 <.0001 <.0001

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Location Estimate  
Year Sixth LA Tech 8.8736 A
Year Sixth      
Year Sixth Hill Farm 5.9524 B

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Third 1 161.5 31.59 <.0001

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Location _Location Estimate Standard Error DF t Value Pr > |t| Adj P
Year Third Hill Farm LA Tech -0.7158 0.1274 161.5 -5.62 <.0001 <.0001

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Location Estimate  
Year Third LA Tech 3.3729 A
Year Third      
Year Third Hill Farm 2.6572 B

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Fifth 5 161.5 17.82 <.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Family _Family Estimate Standard Error DF t Value Pr > |t| Adj P
Year Fifth AGH2 AGH25 -0.5088 0.2206 161.5 -2.31 0.0223 0.1952
Year Fifth AGH2 AGHS1 -1.4087 0.2206 161.5 -6.39 <.0001 <.0001
Year Fifth AGH2 AGHS2 -1.7111 0.2206 161.5 -7.76 <.0001 <.0001
Year Fifth AGH2 AGHS3 -0.7408 0.2206 161.5 -3.36 0.0010 0.0115
Year Fifth AGH2 AGHS4 -1.4383 0.2206 161.5 -6.52 <.0001 <.0001
Year Fifth AGH25 AGHS1 -0.8999 0.2206 161.5 -4.08 <.0001 0.0009
Year Fifth AGH25 AGHS2 -1.2023 0.2206 161.5 -5.45 <.0001 <.0001
Year Fifth AGH25 AGHS3 -0.2320 0.2206 161.5 -1.05 0.2945 0.8997
Year Fifth AGH25 AGHS4 -0.9295 0.2206 161.5 -4.21 <.0001 0.0005
Year Fifth AGHS1 AGHS2 -0.3024 0.2206 161.5 -1.37 0.1723 0.7443
Year Fifth AGHS1 AGHS3 0.6679 0.2206 161.5 3.03 0.0029 0.0321
Year Fifth AGHS1 AGHS4 -0.02959 0.2206 161.5 -0.13 0.8934 1.0000
Year Fifth AGHS2 AGHS3 0.9704 0.2206 161.5 4.40 <.0001 0.0002
Year Fifth AGHS2 AGHS4 0.2728 0.2206 161.5 1.24 0.2179 0.8183
Year Fifth AGHS3 AGHS4 -0.6975 0.2206 161.5 -3.16 0.0019 0.0215

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Family Estimate  
Year Fifth AGHS2 7.1120   A
Year Fifth       A
Year Fifth AGHS4 6.8392   A
Year Fifth       A
Year Fifth AGHS1 6.8096   A
Year Fifth        
Year Fifth AGHS3 6.1416   B
Year Fifth       B
Year Fifth AGH25 5.9097 C B
Year Fifth     C  
Year Fifth AGH2 5.4008 C  

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year First 5 161.5 0.72 0.6060

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Family _Family Estimate Standard Error DF t Value Pr > |t| Adj P
Year First AGH2 AGH25 -0.00289 0.2206 161.5 -0.01 0.9896 1.0000
Year First AGH2 AGHS1 -0.2847 0.2206 161.5 -1.29 0.1986 0.7899
Year First AGH2 AGHS2 -0.2683 0.2206 161.5 -1.22 0.2256 0.8286
Year First AGH2 AGHS3 -0.1805 0.2206 161.5 -0.82 0.4143 0.9639
Year First AGH2 AGHS4 -0.2619 0.2206 161.5 -1.19 0.2368 0.8426
Year First AGH25 AGHS1 -0.2818 0.2206 161.5 -1.28 0.2032 0.7970
Year First AGH25 AGHS2 -0.2654 0.2206 161.5 -1.20 0.2306 0.8350
Year First AGH25 AGHS3 -0.1777 0.2206 161.5 -0.81 0.4218 0.9664
Year First AGH25 AGHS4 -0.2590 0.2206 161.5 -1.17 0.2420 0.8488
Year First AGHS1 AGHS2 0.01642 0.2206 161.5 0.07 0.9408 1.0000
Year First AGHS1 AGHS3 0.1042 0.2206 161.5 0.47 0.6374 0.9971
Year First AGHS1 AGHS4 0.02282 0.2206 161.5 0.10 0.9177 1.0000
Year First AGHS2 AGHS3 0.08776 0.2206 161.5 0.40 0.6913 0.9987
Year First AGHS2 AGHS4 0.006401 0.2206 161.5 0.03 0.9769 1.0000
Year First AGHS3 AGHS4 -0.08136 0.2206 161.5 -0.37 0.7127 0.9991

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Family Estimate  
Year First AGHS1 0.9817 A
Year First     A
Year First AGHS2 0.9653 A
Year First     A
Year First AGHS4 0.9589 A
Year First     A
Year First AGHS3 0.8775 A
Year First     A
Year First AGH25 0.6998 A
Year First     A
Year First AGH2 0.6970 A

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Fourth 5 161.5 20.56 <.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Family _Family Estimate Standard Error DF t Value Pr > |t| Adj P
Year Fourth AGH2 AGH25 -0.3428 0.2206 161.5 -1.55 0.1221 0.6294
Year Fourth AGH2 AGHS1 -1.5691 0.2206 161.5 -7.11 <.0001 <.0001
Year Fourth AGH2 AGHS2 -1.7651 0.2206 161.5 -8.00 <.0001 <.0001
Year Fourth AGH2 AGHS3 -0.9453 0.2206 161.5 -4.29 <.0001 0.0004
Year Fourth AGH2 AGHS4 -1.3959 0.2206 161.5 -6.33 <.0001 <.0001
Year Fourth AGH25 AGHS1 -1.2262 0.2206 161.5 -5.56 <.0001 <.0001
Year Fourth AGH25 AGHS2 -1.4223 0.2206 161.5 -6.45 <.0001 <.0001
Year Fourth AGH25 AGHS3 -0.6025 0.2206 161.5 -2.73 0.0070 0.0727
Year Fourth AGH25 AGHS4 -1.0531 0.2206 161.5 -4.77 <.0001 <.0001
Year Fourth AGHS1 AGHS2 -0.1960 0.2206 161.5 -0.89 0.3754 0.9490
Year Fourth AGHS1 AGHS3 0.6237 0.2206 161.5 2.83 0.0053 0.0564
Year Fourth AGHS1 AGHS4 0.1731 0.2206 161.5 0.78 0.4337 0.9699
Year Fourth AGHS2 AGHS3 0.8198 0.2206 161.5 3.72 0.0003 0.0034
Year Fourth AGHS2 AGHS4 0.3692 0.2206 161.5 1.67 0.0961 0.5505
Year Fourth AGHS3 AGHS4 -0.4506 0.2206 161.5 -2.04 0.0427 0.3212

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Family Estimate  
Year Fourth AGHS2 5.5407   A
Year Fourth       A
Year Fourth AGHS1 5.3447 B A
Year Fourth     B A
Year Fourth AGHS4 5.1715 B A
Year Fourth     B  
Year Fourth AGHS3 4.7209 B C
Year Fourth       C
Year Fourth AGH25 4.1184 D C
Year Fourth     D  
Year Fourth AGH2 3.7756 D  

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Initial 5 161.5 0.14 0.9815

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Family _Family Estimate Standard Error DF t Value Pr > |t| Adj P
Year Initial AGH2 AGH25 0.02700 0.2206 161.5 0.12 0.9027 1.0000
Year Initial AGH2 AGHS1 -0.09858 0.2206 161.5 -0.45 0.6555 0.9977
Year Initial AGH2 AGHS2 -0.02001 0.2206 161.5 -0.09 0.9278 1.0000
Year Initial AGH2 AGHS3 -0.09166 0.2206 161.5 -0.42 0.6783 0.9984
Year Initial AGH2 AGHS4 -0.1142 0.2206 161.5 -0.52 0.6053 0.9955
Year Initial AGH25 AGHS1 -0.1256 0.2206 161.5 -0.57 0.5699 0.9929
Year Initial AGH25 AGHS2 -0.04701 0.2206 161.5 -0.21 0.8315 0.9999
Year Initial AGH25 AGHS3 -0.1187 0.2206 161.5 -0.54 0.5914 0.9946
Year Initial AGH25 AGHS4 -0.1412 0.2206 161.5 -0.64 0.5230 0.9878
Year Initial AGHS1 AGHS2 0.07857 0.2206 161.5 0.36 0.7221 0.9992
Year Initial AGHS1 AGHS3 0.006927 0.2206 161.5 0.03 0.9750 1.0000
Year Initial AGHS1 AGHS4 -0.01562 0.2206 161.5 -0.07 0.9436 1.0000
Year Initial AGHS2 AGHS3 -0.07165 0.2206 161.5 -0.32 0.7457 0.9995
Year Initial AGHS2 AGHS4 -0.09420 0.2206 161.5 -0.43 0.6699 0.9982
Year Initial AGHS3 AGHS4 -0.02255 0.2206 161.5 -0.10 0.9187 1.0000

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Family Estimate  
Year Initial AGHS4 0.5520 A
Year Initial     A
Year Initial AGHS1 0.5363 A
Year Initial     A
Year Initial AGHS3 0.5294 A
Year Initial     A
Year Initial AGHS2 0.4578 A
Year Initial     A
Year Initial AGH2 0.4378 A
Year Initial     A
Year Initial AGH25 0.4108 A

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Second 5 161.5 5.65 <.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Family _Family Estimate Standard Error DF t Value Pr > |t| Adj P
Year Second AGH2 AGH25 -0.09287 0.2206 161.5 -0.42 0.6743 0.9983
Year Second AGH2 AGHS1 -0.7712 0.2206 161.5 -3.50 0.0006 0.0073
Year Second AGH2 AGHS2 -0.9229 0.2206 161.5 -4.18 <.0001 0.0006
Year Second AGH2 AGHS3 -0.5018 0.2206 161.5 -2.28 0.0242 0.2082
Year Second AGH2 AGHS4 -0.6426 0.2206 161.5 -2.91 0.0041 0.0446
Year Second AGH25 AGHS1 -0.6784 0.2206 161.5 -3.08 0.0025 0.0279
Year Second AGH25 AGHS2 -0.8300 0.2206 161.5 -3.76 0.0002 0.0028
Year Second AGH25 AGHS3 -0.4090 0.2206 161.5 -1.85 0.0656 0.4332
Year Second AGH25 AGHS4 -0.5497 0.2206 161.5 -2.49 0.0137 0.1303
Year Second AGHS1 AGHS2 -0.1517 0.2206 161.5 -0.69 0.4927 0.9832
Year Second AGHS1 AGHS3 0.2694 0.2206 161.5 1.22 0.2237 0.8261
Year Second AGHS1 AGHS4 0.1286 0.2206 161.5 0.58 0.5606 0.9921
Year Second AGHS2 AGHS3 0.4211 0.2206 161.5 1.91 0.0580 0.3991
Year Second AGHS2 AGHS4 0.2803 0.2206 161.5 1.27 0.2056 0.8007
Year Second AGHS3 AGHS4 -0.1408 0.2206 161.5 -0.64 0.5243 0.9880

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Family Estimate  
Year Second AGHS2 2.5035   A  
Year Second       A  
Year Second AGHS1 2.3519   A  
Year Second       A  
Year Second AGHS4 2.2232 B A  
Year Second     B A  
Year Second AGHS3 2.0825 B A C
Year Second     B   C
Year Second AGH25 1.6735 B   C
Year Second         C
Year Second AGH2 1.5806     C

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Seventh 5 161.5 15.75 <.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Family _Family Estimate Standard Error DF t Value Pr > |t| Adj P
Year Seventh AGH2 AGH25 0.1012 0.2206 161.5 0.46 0.6469 0.9974
Year Seventh AGH2 AGHS1 -1.2296 0.2206 161.5 -5.57 <.0001 <.0001
Year Seventh AGH2 AGHS2 -1.2242 0.2206 161.5 -5.55 <.0001 <.0001
Year Seventh AGH2 AGHS3 -0.5381 0.2206 161.5 -2.44 0.0158 0.1468
Year Seventh AGH2 AGHS4 -1.1547 0.2206 161.5 -5.23 <.0001 <.0001
Year Seventh AGH25 AGHS1 -1.3308 0.2206 161.5 -6.03 <.0001 <.0001
Year Seventh AGH25 AGHS2 -1.3254 0.2206 161.5 -6.01 <.0001 <.0001
Year Seventh AGH25 AGHS3 -0.6393 0.2206 161.5 -2.90 0.0043 0.0465
Year Seventh AGH25 AGHS4 -1.2559 0.2206 161.5 -5.69 <.0001 <.0001
Year Seventh AGHS1 AGHS2 0.005355 0.2206 161.5 0.02 0.9807 1.0000
Year Seventh AGHS1 AGHS3 0.6915 0.2206 161.5 3.13 0.0020 0.0234
Year Seventh AGHS1 AGHS4 0.07490 0.2206 161.5 0.34 0.7346 0.9994
Year Seventh AGHS2 AGHS3 0.6861 0.2206 161.5 3.11 0.0022 0.0251
Year Seventh AGHS2 AGHS4 0.06954 0.2206 161.5 0.32 0.7530 0.9996
Year Seventh AGHS3 AGHS4 -0.6166 0.2206 161.5 -2.80 0.0058 0.0615

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Family Estimate  
Year Seventh AGHS1 9.6536   A
Year Seventh       A
Year Seventh AGHS2 9.6482   A
Year Seventh       A
Year Seventh AGHS4 9.5787 B A
Year Seventh     B  
Year Seventh AGHS3 8.9621 B C
Year Seventh       C
Year Seventh AGH2 8.4240 D C
Year Seventh     D  
Year Seventh AGH25 8.3228 D  

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Sixth 5 161.5 27.96 <.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Family _Family Estimate Standard Error DF t Value Pr > |t| Adj P
Year Sixth AGH2 AGH25 -0.3353 0.2206 161.5 -1.52 0.1304 0.6515
Year Sixth AGH2 AGHS1 -1.8840 0.2206 161.5 -8.54 <.0001 <.0001
Year Sixth AGH2 AGHS2 -1.9521 0.2206 161.5 -8.85 <.0001 <.0001
Year Sixth AGH2 AGHS3 -1.2217 0.2206 161.5 -5.54 <.0001 <.0001
Year Sixth AGH2 AGHS4 -1.6303 0.2206 161.5 -7.39 <.0001 <.0001
Year Sixth AGH25 AGHS1 -1.5486 0.2206 161.5 -7.02 <.0001 <.0001
Year Sixth AGH25 AGHS2 -1.6167 0.2206 161.5 -7.33 <.0001 <.0001
Year Sixth AGH25 AGHS3 -0.8864 0.2206 161.5 -4.02 <.0001 0.0011
Year Sixth AGH25 AGHS4 -1.2950 0.2206 161.5 -5.87 <.0001 <.0001
Year Sixth AGHS1 AGHS2 -0.06811 0.2206 161.5 -0.31 0.7579 0.9996
Year Sixth AGHS1 AGHS3 0.6623 0.2206 161.5 3.00 0.0031 0.0346
Year Sixth AGHS1 AGHS4 0.2536 0.2206 161.5 1.15 0.2519 0.8599
Year Sixth AGHS2 AGHS3 0.7304 0.2206 161.5 3.31 0.0011 0.0135
Year Sixth AGHS2 AGHS4 0.3217 0.2206 161.5 1.46 0.1466 0.6909
Year Sixth AGHS3 AGHS4 -0.4086 0.2206 161.5 -1.85 0.0658 0.4341

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Family Estimate  
Year Sixth AGHS2 8.1945   A
Year Sixth       A
Year Sixth AGHS1 8.1264   A
Year Sixth       A
Year Sixth AGHS4 7.8728 B A
Year Sixth     B  
Year Sixth AGHS3 7.4642 B  
Year Sixth        
Year Sixth AGH25 6.5778   C
Year Sixth       C
Year Sixth AGH2 6.2425   C

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Third 5 161.5 8.12 <.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Family _Family Estimate Standard Error DF t Value Pr > |t| Adj P
Year Third AGH2 AGH25 -0.1341 0.2206 161.5 -0.61 0.5442 0.9904
Year Third AGH2 AGHS1 -1.0244 0.2206 161.5 -4.64 <.0001 <.0001
Year Third AGH2 AGHS2 -1.0013 0.2206 161.5 -4.54 <.0001 0.0001
Year Third AGH2 AGHS3 -0.3330 0.2206 161.5 -1.51 0.1331 0.6585
Year Third AGH2 AGHS4 -0.7438 0.2206 161.5 -3.37 0.0009 0.0110
Year Third AGH25 AGHS1 -0.8903 0.2206 161.5 -4.04 <.0001 0.0010
Year Third AGH25 AGHS2 -0.8673 0.2206 161.5 -3.93 0.0001 0.0015
Year Third AGH25 AGHS3 -0.1989 0.2206 161.5 -0.90 0.3686 0.9458
Year Third AGH25 AGHS4 -0.6097 0.2206 161.5 -2.76 0.0064 0.0667
Year Third AGHS1 AGHS2 0.02302 0.2206 161.5 0.10 0.9170 1.0000
Year Third AGHS1 AGHS3 0.6914 0.2206 161.5 3.13 0.0020 0.0234
Year Third AGHS1 AGHS4 0.2806 0.2206 161.5 1.27 0.2052 0.8000
Year Third AGHS2 AGHS3 0.6684 0.2206 161.5 3.03 0.0028 0.0319
Year Third AGHS2 AGHS4 0.2576 0.2206 161.5 1.17 0.2447 0.8518
Year Third AGHS3 AGHS4 -0.4108 0.2206 161.5 -1.86 0.0643 0.4278

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Family Estimate  
Year Third AGHS1 3.5000   A
Year Third       A
Year Third AGHS2 3.4770   A
Year Third       A
Year Third AGHS4 3.2194 B A
Year Third     B  
Year Third AGHS3 2.8086 B C
Year Third     B C
Year Third AGH25 2.6097 B C
Year Third       C
Year Third AGH2 2.4756   C

LS-Means Results

Family LS-Means

Family Least Squares Means
Family Estimate Standard
Error
DF t Value Pr > |t|
AGH2 3.6292 0.1008 36 36.01 <.0001
AGH25 3.7903 0.1008 36 37.60 <.0001
AGHS1 4.6630 0.1008 36 46.26 <.0001
AGHS2 4.7374 0.1008 36 47.00 <.0001
AGHS3 4.1983 0.1008 36 41.65 <.0001
AGHS4 4.5520 0.1008 36 45.16 <.0001

Family Diffs

Differences of Family Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Family _Family Estimate Standard Error DF t Value Pr > |t| Adj P
AGH2 AGH25 -0.1611 0.1425 36 -1.13 0.2660 0.8656
AGH2 AGHS1 -1.0338 0.1425 36 -7.25 <.0001 <.0001
AGH2 AGHS2 -1.1081 0.1425 36 -7.77 <.0001 <.0001
AGH2 AGHS3 -0.5691 0.1425 36 -3.99 0.0003 0.0039
AGH2 AGHS4 -0.9227 0.1425 36 -6.47 <.0001 <.0001
AGH25 AGHS1 -0.8727 0.1425 36 -6.12 <.0001 <.0001
AGH25 AGHS2 -0.9471 0.1425 36 -6.64 <.0001 <.0001
AGH25 AGHS3 -0.4080 0.1425 36 -2.86 0.0070 0.0700
AGH25 AGHS4 -0.7616 0.1425 36 -5.34 <.0001 <.0001
AGHS1 AGHS2 -0.07436 0.1425 36 -0.52 0.6051 0.9950
AGHS1 AGHS3 0.4647 0.1425 36 3.26 0.0024 0.0272
AGHS1 AGHS4 0.1111 0.1425 36 0.78 0.4410 0.9693
AGHS2 AGHS3 0.5390 0.1425 36 3.78 0.0006 0.0069
AGHS2 AGHS4 0.1854 0.1425 36 1.30 0.2016 0.7827
AGHS3 AGHS4 -0.3536 0.1425 36 -2.48 0.0179 0.1569

Family Lines

Tukey Grouping for Family Least Squares Means (Alpha=0.05)
LS-means with the same letter are not significantly different.
Family Estimate  
AGHS2 4.7374   A
      A
AGHS1 4.6630   A
      A
AGHS4 4.5520 B A
    B  
AGHS3 4.1983 B C
      C
AGH25 3.7903 D C
    D  
AGH2 3.6292 D  

SAS Output Crown

htmltools::includeHTML("https://raw.githubusercontent.com/azd169/SAS_outputs/main/hane-results-crown.html")
## Warning: `includeHTML()` was provided a `path` that appears to be a complete HTML document.
## ✖ Path: https://raw.githubusercontent.com/azd169/SAS_outputs/main/hane-results-crown.html
## ℹ Use `tags$iframe()` to include an HTML document. You can either ensure `path` is accessible in your app or document (see e.g. `shiny::addResourcePath()`) and pass the relative path to the `src` argument. Or you can read the contents of `path` and pass the contents to `srcdoc`.
Results: hane.sas

Results: hane.sas

The GLIMMIX Procedure

The Glimmix Procedure

Model Information

Model Information
Data Set WORK.CROWN
Response Variable Crown
Response Distribution Gaussian
Link Function Identity
Variance Function Default
Variance Matrix Blocked By GLB
Estimation Technique Restricted Maximum Likelihood
Degrees of Freedom Method Kenward-Roger
Fixed Effects SE Adjustment Kenward-Roger

Class Level Information

Class Level Information
Class Levels Values
Year 3 Fifth Seventh Sixth
Location 2 Hill Farm LA Tech
Block 4 1 2 3 4
Family 6 AGH2 AGH25 AGHS1 AGHS2 AGHS3 AGHS4
GLB 48 AGH25_Hill Farm_1 AGH25_Hill Farm_2 AGH25_Hill Farm_3 AGH25_Hill Farm_4 AGH25_LA Tech_1 AGH25_LA Tech_2 AGH25_LA Tech_3 AGH25_LA Tech_4 AGH2_Hill Farm_1 AGH2_Hill Farm_2 AGH2_Hill Farm_3 AGH2_Hill Farm_4 AGH2_LA Tech_1 AGH2_LA Tech_2 AGH2_LA Tech_3 AGH2_LA Tech_4 AGHS1_Hill Farm_1 AGHS1_Hill Farm_2 AGHS1_Hill Farm_3 AGHS1_Hill Farm_4 AGHS1_LA Tech_1 AGHS1_LA Tech_2 AGHS1_LA Tech_3 AGHS1_LA Tech_4 AGHS2_Hill Farm_1 AGHS2_Hill Farm_2 AGHS2_Hill Farm_3 AGHS2_Hill Farm_4 AGHS2_LA Tech_1 AGHS2_LA Tech_2 AGHS2_LA Tech_3 AGHS2_LA Tech_4 AGHS3_Hill Farm_1 AGHS3_Hill Farm_2 AGHS3_Hill Farm_3 AGHS3_Hill Farm_4 AGHS3_LA Tech_1 AGHS3_LA Tech_2 AGHS3_LA Tech_3 AGHS3_LA Tech_4 AGHS4_Hill Farm_1 AGHS4_Hill Farm_2 AGHS4_Hill Farm_3 AGHS4_Hill Farm_4 AGHS4_LA Tech_1 AGHS4_LA Tech_2 AGHS4_LA Tech_3 AGHS4_LA Tech_4

Number of Observations

Number of Observations Read 144
Number of Observations Used 144

Dimensions

Dimensions
G-side Cov. Parameters 2
R-side Cov. Parameters 1
Columns in X 84
Columns in Z per Subject 1
Subjects (Blocks in V) 48
Max Obs per Subject 3

Optimization Information

Optimization Information
Optimization Technique Dual Quasi-Newton
Parameters in Optimization 2
Lower Boundaries 2
Upper Boundaries 1
Fixed Effects Profiled
Residual Variance Profiled
Starting From Data

Iteration History

Iteration History
Iteration Restarts Evaluations Objective
Function
Change Max
Gradient
0 0 4 327.11392315 . 7.57E-13

Convergence Status

Convergence criterion (ABSGCONV=0.00001) satisfied.

Hessian Matrix

Hessian Matrix
Cov Parm Subject Row Col1 Col2
Variance GLB 1 14.9172  
AR(1) GLB 2    

Fit Statistics

Fit Statistics
-2 Res Log Likelihood 327.11
AIC (smaller is better) 333.11
AICC (smaller is better) 333.34
BIC (smaller is better) 338.73
CAIC (smaller is better) 341.73
HQIC (smaller is better) 335.24
Generalized Chi-Square 52.75
Gener. Chi-Square / DF 0.49

Covariance Parameter Estimates

Covariance Parameter Estimates
Cov Parm Subject Estimate Standard
Error
Variance GLB 0.4567 0.1485
AR(1) GLB 0 .
Residual   0.4884 0.08140

Solutions for Fixed Effects

Solutions for Fixed Effects
Effect Year Location Family Estimate Standard
Error
DF t Value Pr > |t|
Intercept       4.6956 0.4861 73.62 9.66 <.0001
Year Fifth     -3.0771 0.4942 72 -6.23 <.0001
Year Seventh     1.0079 0.4942 72 2.04 0.0451
Year Sixth     0 . . . .
Location   Hill Farm   -2.5922 0.6874 73.62 -3.77 0.0003
Location   LA Tech   0 . . . .
Year*Location Fifth Hill Farm   1.4553 0.6989 72 2.08 0.0409
Year*Location Fifth LA Tech   0 . . . .
Year*Location Seventh Hill Farm   -0.1176 0.6989 72 -0.17 0.8668
Year*Location Seventh LA Tech   0 . . . .
Year*Location Sixth Hill Farm   0 . . . .
Year*Location Sixth LA Tech   0 . . . .
Family     AGH2 -3.0492 0.6874 73.62 -4.44 <.0001
Family     AGH25 -2.3896 0.6874 73.62 -3.48 0.0009
Family     AGHS1 0.5941 0.6874 73.62 0.86 0.3902
Family     AGHS2 0.8519 0.6874 73.62 1.24 0.2192
Family     AGHS3 -0.08741 0.6874 73.62 -0.13 0.8992
Family     AGHS4 0 . . . .
Year*Family Fifth   AGH2 2.6497 0.6989 72 3.79 0.0003
Year*Family Fifth   AGH25 1.7974 0.6989 72 2.57 0.0122
Year*Family Fifth   AGHS1 -0.1967 0.6989 72 -0.28 0.7791
Year*Family Fifth   AGHS2 0.2015 0.6989 72 0.29 0.7739
Year*Family Fifth   AGHS3 0.3601 0.6989 72 0.52 0.6079
Year*Family Fifth   AGHS4 0 . . . .
Year*Family Seventh   AGH2 0.2702 0.6989 72 0.39 0.7002
Year*Family Seventh   AGH25 -0.4015 0.6989 72 -0.57 0.5674
Year*Family Seventh   AGHS1 -0.2398 0.6989 72 -0.34 0.7325
Year*Family Seventh   AGHS2 -0.3092 0.6989 72 -0.44 0.6595
Year*Family Seventh   AGHS3 -0.6890 0.6989 72 -0.99 0.3275
Year*Family Seventh   AGHS4 0 . . . .
Year*Family Sixth   AGH2 0 . . . .
Year*Family Sixth   AGH25 0 . . . .
Year*Family Sixth   AGHS1 0 . . . .
Year*Family Sixth   AGHS2 0 . . . .
Year*Family Sixth   AGHS3 0 . . . .
Year*Family Sixth   AGHS4 0 . . . .
Location*Family   Hill Farm AGH2 1.8355 0.9722 73.62 1.89 0.0630
Location*Family   Hill Farm AGH25 1.4387 0.9722 73.62 1.48 0.1432
Location*Family   Hill Farm AGHS1 0.1333 0.9722 73.62 0.14 0.8913
Location*Family   Hill Farm AGHS2 0.4140 0.9722 73.62 0.43 0.6714
Location*Family   Hill Farm AGHS3 0.1992 0.9722 73.62 0.20 0.8382
Location*Family   Hill Farm AGHS4 0 . . . .
Location*Family   LA Tech AGH2 0 . . . .
Location*Family   LA Tech AGH25 0 . . . .
Location*Family   LA Tech AGHS1 0 . . . .
Location*Family   LA Tech AGHS2 0 . . . .
Location*Family   LA Tech AGHS3 0 . . . .
Location*Family   LA Tech AGHS4 0 . . . .
Year*Location*Family Fifth Hill Farm AGH2 -1.4704 0.9883 72 -1.49 0.1412
Year*Location*Family Fifth Hill Farm AGH25 -0.8170 0.9883 72 -0.83 0.4112
Year*Location*Family Fifth Hill Farm AGHS1 -0.1838 0.9883 72 -0.19 0.8530
Year*Location*Family Fifth Hill Farm AGHS2 -1.0488 0.9883 72 -1.06 0.2921
Year*Location*Family Fifth Hill Farm AGHS3 -0.4048 0.9883 72 -0.41 0.6833
Year*Location*Family Fifth Hill Farm AGHS4 0 . . . .
Year*Location*Family Fifth LA Tech AGH2 0 . . . .
Year*Location*Family Fifth LA Tech AGH25 0 . . . .
Year*Location*Family Fifth LA Tech AGHS1 0 . . . .
Year*Location*Family Fifth LA Tech AGHS2 0 . . . .
Year*Location*Family Fifth LA Tech AGHS3 0 . . . .
Year*Location*Family Fifth LA Tech AGHS4 0 . . . .
Year*Location*Family Seventh Hill Farm AGH2 -0.5397 0.9883 72 -0.55 0.5867
Year*Location*Family Seventh Hill Farm AGH25 -0.4144 0.9883 72 -0.42 0.6762
Year*Location*Family Seventh Hill Farm AGHS1 0.1666 0.9883 72 0.17 0.8666
Year*Location*Family Seventh Hill Farm AGHS2 0.3237 0.9883 72 0.33 0.7443
Year*Location*Family Seventh Hill Farm AGHS3 0.5593 0.9883 72 0.57 0.5732
Year*Location*Family Seventh Hill Farm AGHS4 0 . . . .
Year*Location*Family Seventh LA Tech AGH2 0 . . . .
Year*Location*Family Seventh LA Tech AGH25 0 . . . .
Year*Location*Family Seventh LA Tech AGHS1 0 . . . .
Year*Location*Family Seventh LA Tech AGHS2 0 . . . .
Year*Location*Family Seventh LA Tech AGHS3 0 . . . .
Year*Location*Family Seventh LA Tech AGHS4 0 . . . .
Year*Location*Family Sixth Hill Farm AGH2 0 . . . .
Year*Location*Family Sixth Hill Farm AGH25 0 . . . .
Year*Location*Family Sixth Hill Farm AGHS1 0 . . . .
Year*Location*Family Sixth Hill Farm AGHS2 0 . . . .
Year*Location*Family Sixth Hill Farm AGHS3 0 . . . .
Year*Location*Family Sixth Hill Farm AGHS4 0 . . . .
Year*Location*Family Sixth LA Tech AGH2 0 . . . .
Year*Location*Family Sixth LA Tech AGH25 0 . . . .
Year*Location*Family Sixth LA Tech AGHS1 0 . . . .
Year*Location*Family Sixth LA Tech AGHS2 0 . . . .
Year*Location*Family Sixth LA Tech AGHS3 0 . . . .
Year*Location*Family Sixth LA Tech AGHS4 0 . . . .

Type III Tests of Fixed Effects

Type III Tests of Fixed Effects
Effect Num DF Den DF F Value Pr > F
Year 2 72 177.27 <.0001
Location 1 36 55.25 <.0001
Year*Location 2 72 6.01 0.0039
Family 5 36 12.64 <.0001
Year*Family 10 72 5.39 <.0001
Location*Family 5 36 0.82 0.5439
Year*Location*Family 10 72 0.51 0.8776

Solution for Random Effects

Solution for Random Effects
Effect Subject Estimate Std Err Pred DF t Value Pr > |t| Alpha Lower Upper
Intercept GLB AGH25_Hill Farm_1 -0.1422 0.4578 42.76 -0.31 0.7575 0.05 -1.0655 0.7811
Intercept GLB AGH25_Hill Farm_2 0.2869 0.4578 42.76 0.63 0.5341 0.05 -0.6364 1.2102
Intercept GLB AGH25_Hill Farm_3 -0.2848 0.4578 42.76 -0.62 0.5371 0.05 -1.2081 0.6385
Intercept GLB AGH25_Hill Farm_4 0.1401 0.4578 42.76 0.31 0.7611 0.05 -0.7832 1.0634
Intercept GLB AGH25_LA Tech_1 -0.6896 0.4578 42.76 -1.51 0.1393 0.05 -1.6129 0.2337
Intercept GLB AGH25_LA Tech_2 0.1718 0.4578 42.76 0.38 0.7093 0.05 -0.7515 1.0951
Intercept GLB AGH25_LA Tech_3 0.2293 0.4578 42.76 0.50 0.6190 0.05 -0.6940 1.1526
Intercept GLB AGH25_LA Tech_4 0.2885 0.4578 42.76 0.63 0.5319 0.05 -0.6348 1.2118
Intercept GLB AGH2_Hill Farm_1 -0.3502 0.4578 42.76 -0.77 0.4484 0.05 -1.2735 0.5731
Intercept GLB AGH2_Hill Farm_2 0.2129 0.4578 42.76 0.47 0.6442 0.05 -0.7104 1.1362
Intercept GLB AGH2_Hill Farm_3 0.2285 0.4578 42.76 0.50 0.6201 0.05 -0.6948 1.1518
Intercept GLB AGH2_Hill Farm_4 -0.09127 0.4578 42.76 -0.20 0.8429 0.05 -1.0146 0.8320
Intercept GLB AGH2_LA Tech_1 -0.5441 0.4578 42.76 -1.19 0.2412 0.05 -1.4673 0.3792
Intercept GLB AGH2_LA Tech_2 -0.2853 0.4578 42.76 -0.62 0.5364 0.05 -1.2086 0.6380
Intercept GLB AGH2_LA Tech_3 0.8694 0.4578 42.76 1.90 0.0643 0.05 -0.05387 1.7927
Intercept GLB AGH2_LA Tech_4 -0.04006 0.4578 42.76 -0.09 0.9307 0.05 -0.9634 0.8832
Intercept GLB AGHS1_Hill Farm_1 -0.4909 0.4578 42.76 -1.07 0.2896 0.05 -1.4142 0.4324
Intercept GLB AGHS1_Hill Farm_2 0.5528 0.4578 42.76 1.21 0.2338 0.05 -0.3705 1.4761
Intercept GLB AGHS1_Hill Farm_3 0.4166 0.4578 42.76 0.91 0.3679 0.05 -0.5067 1.3399
Intercept GLB AGHS1_Hill Farm_4 -0.4785 0.4578 42.76 -1.05 0.3017 0.05 -1.4018 0.4448
Intercept GLB AGHS1_LA Tech_1 -1.1449 0.4578 42.76 -2.50 0.0163 0.05 -2.0682 -0.2216
Intercept GLB AGHS1_LA Tech_2 -0.2571 0.4578 42.76 -0.56 0.5772 0.05 -1.1804 0.6662
Intercept GLB AGHS1_LA Tech_3 0.8501 0.4578 42.76 1.86 0.0702 0.05 -0.07321 1.7734
Intercept GLB AGHS1_LA Tech_4 0.5519 0.4578 42.76 1.21 0.2346 0.05 -0.3714 1.4752
Intercept GLB AGHS2_Hill Farm_1 -0.9490 0.4578 42.76 -2.07 0.0442 0.05 -1.8723 -0.02572
Intercept GLB AGHS2_Hill Farm_2 0.3246 0.4578 42.76 0.71 0.4821 0.05 -0.5987 1.2479
Intercept GLB AGHS2_Hill Farm_3 0.2169 0.4578 42.76 0.47 0.6381 0.05 -0.7064 1.1402
Intercept GLB AGHS2_Hill Farm_4 0.4075 0.4578 42.76 0.89 0.3783 0.05 -0.5158 1.3308
Intercept GLB AGHS2_LA Tech_1 -0.9095 0.4578 42.76 -1.99 0.0534 0.05 -1.8328 0.01377
Intercept GLB AGHS2_LA Tech_2 0.06896 0.4578 42.76 0.15 0.8810 0.05 -0.8543 0.9923
Intercept GLB AGHS2_LA Tech_3 0.7353 0.4578 42.76 1.61 0.1156 0.05 -0.1880 1.6586
Intercept GLB AGHS2_LA Tech_4 0.1053 0.4578 42.76 0.23 0.8192 0.05 -0.8180 1.0286
Intercept GLB AGHS3_Hill Farm_1 -0.4208 0.4578 42.76 -0.92 0.3631 0.05 -1.3441 0.5025
Intercept GLB AGHS3_Hill Farm_2 -0.4823 0.4578 42.76 -1.05 0.2980 0.05 -1.4056 0.4410
Intercept GLB AGHS3_Hill Farm_3 0.6032 0.4578 42.76 1.32 0.1946 0.05 -0.3201 1.5265
Intercept GLB AGHS3_Hill Farm_4 0.2999 0.4578 42.76 0.66 0.5159 0.05 -0.6234 1.2232
Intercept GLB AGHS3_LA Tech_1 -0.6954 0.4578 42.76 -1.52 0.1361 0.05 -1.6187 0.2279
Intercept GLB AGHS3_LA Tech_2 -0.3490 0.4578 42.76 -0.76 0.4500 0.05 -1.2723 0.5743
Intercept GLB AGHS3_LA Tech_3 0.9425 0.4578 42.76 2.06 0.0456 0.05 0.01919 1.8658
Intercept GLB AGHS3_LA Tech_4 0.1019 0.4578 42.76 0.22 0.8249 0.05 -0.8214 1.0252
Intercept GLB AGHS4_Hill Farm_1 -0.7912 0.4578 42.76 -1.73 0.0911 0.05 -1.7145 0.1321
Intercept GLB AGHS4_Hill Farm_2 0.06875 0.4578 42.76 0.15 0.8813 0.05 -0.8545 0.9920
Intercept GLB AGHS4_Hill Farm_3 0.005591 0.4578 42.76 0.01 0.9903 0.05 -0.9177 0.9289
Intercept GLB AGHS4_Hill Farm_4 0.7168 0.4578 42.76 1.57 0.1247 0.05 -0.2065 1.6401
Intercept GLB AGHS4_LA Tech_1 -0.5740 0.4578 42.76 -1.25 0.2167 0.05 -1.4973 0.3493
Intercept GLB AGHS4_LA Tech_2 0.2233 0.4578 42.76 0.49 0.6281 0.05 -0.7000 1.1466
Intercept GLB AGHS4_LA Tech_3 0.2200 0.4578 42.76 0.48 0.6332 0.05 -0.7033 1.1433
Intercept GLB AGHS4_LA Tech_4 0.1306 0.4578 42.76 0.29 0.7768 0.05 -0.7927 1.0539

Slices

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Fifth 5 73.62 1.24 0.2982

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Family _Family Estimate Standard Error DF t Value Pr > |t| Adj P
Year Fifth AGH2 AGH25 0.06438 0.4861 73.62 0.13 0.8950 1.0000
Year Fifth AGH2 AGHS1 -0.5890 0.4861 73.62 -1.21 0.2295 0.8297
Year Fifth AGH2 AGHS2 -0.9529 0.4861 73.62 -1.96 0.0537 0.3750
Year Fifth AGH2 AGHS3 -0.3868 0.4861 73.62 -0.80 0.4287 0.9674
Year Fifth AGH2 AGHS4 -0.2169 0.4861 73.62 -0.45 0.6567 0.9977
Year Fifth AGH25 AGHS1 -0.6534 0.4861 73.62 -1.34 0.1830 0.7594
Year Fifth AGH25 AGHS2 -1.0173 0.4861 73.62 -2.09 0.0398 0.3027
Year Fifth AGH25 AGHS3 -0.4512 0.4861 73.62 -0.93 0.3563 0.9379
Year Fifth AGH25 AGHS4 -0.2813 0.4861 73.62 -0.58 0.5646 0.9921
Year Fifth AGHS1 AGHS2 -0.3639 0.4861 73.62 -0.75 0.4565 0.9750
Year Fifth AGHS1 AGHS3 0.2023 0.4861 73.62 0.42 0.6786 0.9983
Year Fifth AGHS1 AGHS4 0.3721 0.4861 73.62 0.77 0.4464 0.9724
Year Fifth AGHS2 AGHS3 0.5661 0.4861 73.62 1.16 0.2479 0.8520
Year Fifth AGHS2 AGHS4 0.7360 0.4861 73.62 1.51 0.1343 0.6564
Year Fifth AGHS3 AGHS4 0.1699 0.4861 73.62 0.35 0.7277 0.9993

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Family Estimate  
Year Fifth AGHS2 1.7861 A
Year Fifth     A
Year Fifth AGHS1 1.4222 A
Year Fifth     A
Year Fifth AGHS3 1.2199 A
Year Fifth     A
Year Fifth AGHS4 1.0501 A
Year Fifth     A
Year Fifth AGH2 0.8332 A
Year Fifth     A
Year Fifth AGH25 0.7688 A

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Seventh 5 73.62 15.34 <.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Family _Family Estimate Standard Error DF t Value Pr > |t| Adj P
Year Seventh AGH2 AGH25 0.1479 0.4861 73.62 0.30 0.7619 0.9996
Year Seventh AGH2 AGHS1 -2.6354 0.4861 73.62 -5.42 <.0001 <.0001
Year Seventh AGH2 AGHS2 -3.0426 0.4861 73.62 -6.26 <.0001 <.0001
Year Seventh AGH2 AGHS3 -1.7339 0.4861 73.62 -3.57 0.0006 0.0082
Year Seventh AGH2 AGHS4 -2.1311 0.4861 73.62 -4.38 <.0001 0.0005
Year Seventh AGH25 AGHS1 -2.7833 0.4861 73.62 -5.73 <.0001 <.0001
Year Seventh AGH25 AGHS2 -3.1904 0.4861 73.62 -6.56 <.0001 <.0001
Year Seventh AGH25 AGHS3 -1.8818 0.4861 73.62 -3.87 0.0002 0.0031
Year Seventh AGH25 AGHS4 -2.2789 0.4861 73.62 -4.69 <.0001 0.0002
Year Seventh AGHS1 AGHS2 -0.4072 0.4861 73.62 -0.84 0.4049 0.9594
Year Seventh AGHS1 AGHS3 0.9015 0.4861 73.62 1.85 0.0677 0.4382
Year Seventh AGHS1 AGHS4 0.5043 0.4861 73.62 1.04 0.3029 0.9036
Year Seventh AGHS2 AGHS3 1.3087 0.4861 73.62 2.69 0.0088 0.0893
Year Seventh AGHS2 AGHS4 0.9115 0.4861 73.62 1.88 0.0647 0.4256
Year Seventh AGHS3 AGHS4 -0.3972 0.4861 73.62 -0.82 0.4165 0.9635

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Family Estimate  
Year Seventh AGHS2 5.2601 A
Year Seventh     A
Year Seventh AGHS1 4.8530 A
Year Seventh     A
Year Seventh AGHS4 4.3486 A
Year Seventh     A
Year Seventh AGHS3 3.9514 A
Year Seventh      
Year Seventh AGH2 2.2175 B
Year Seventh     B
Year Seventh AGH25 2.0697 B

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice Num DF Den DF F Value Pr > F
Year Sixth 5 73.62 13.84 <.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means
Adjustment for Multiple Comparisons: Tukey
Slice Family _Family Estimate Standard Error DF t Value Pr > |t| Adj P
Year Sixth AGH2 AGH25 -0.4612 0.4861 73.62 -0.95 0.3459 0.9322
Year Sixth AGH2 AGHS1 -2.7922 0.4861 73.62 -5.74 <.0001 <.0001
Year Sixth AGH2 AGHS2 -3.1903 0.4861 73.62 -6.56 <.0001 <.0001
Year Sixth AGH2 AGHS3 -2.1436 0.4861 73.62 -4.41 <.0001 0.0005
Year Sixth AGH2 AGHS4 -2.1314 0.4861 73.62 -4.38 <.0001 0.0005
Year Sixth AGH25 AGHS1 -2.3311 0.4861 73.62 -4.80 <.0001 0.0001
Year Sixth AGH25 AGHS2 -2.7291 0.4861 73.62 -5.61 <.0001 <.0001
Year Sixth AGH25 AGHS3 -1.6824 0.4861 73.62 -3.46 0.0009 0.0113
Year Sixth AGH25 AGHS4 -1.6703 0.4861 73.62 -3.44 0.0010 0.0122
Year Sixth AGHS1 AGHS2 -0.3981 0.4861 73.62 -0.82 0.4155 0.9631
Year Sixth AGHS1 AGHS3 0.6486 0.4861 73.62 1.33 0.1862 0.7650
Year Sixth AGHS1 AGHS4 0.6608 0.4861 73.62 1.36 0.1782 0.7507
Year Sixth AGHS2 AGHS3 1.0467 0.4861 73.62 2.15 0.0346 0.2725
Year Sixth AGHS2 AGHS4 1.0589 0.4861 73.62 2.18 0.0326 0.2606
Year Sixth AGHS3 AGHS4 0.01218 0.4861 73.62 0.03 0.9801 1.0000

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice Family Estimate  
Year Sixth AGHS2 4.4584 A
Year Sixth     A
Year Sixth AGHS1 4.0603 A
Year Sixth     A
Year Sixth AGHS3 3.4117 A
Year Sixth     A
Year Sixth AGHS4 3.3995 A
Year Sixth      
Year Sixth AGH25 1.7292 B
Year Sixth     B
Year Sixth AGH2 1.2681 B

Survival rates

Sur1 <- ggplot(hane_survival_year,
               aes(x = Survival, y = Rate, fill = Survival, Color = Survival)) +
  stat_summary(fun = mean, geom = "bar", position = "dodge") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 100), breaks = c(0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("green4", "red3"), 0.5)) +
  labs(title = "Survival Rate by Growing Season") +
  xlab("") +
  ylab("Survival rate (%)") +
  facet_wrap( ~ Year) +
  guides(fill = "none") +
  theme_classic()

plotly::ggplotly(Sur1)
Sur2 <- ggplot(hane_survival_year_family,
               aes(x = Family, y = Rate, fill = Survival, color = Survival)) +
  stat_summary(fun = mean, geom = "bar", position = "dodge") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 100), breaks = c(0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("green4", "red3"), 0.5)) +
  labs(title = "Survival Rate by Growing Season and Family") +
  xlab("") +
  ylab("Survival rate (%)") +
  facet_wrap( ~ Year) +
  guides(fill = "none") +
  theme_classic()

plotly::ggplotly(Sur2)
Sur3 <- ggplot(hane_survival_year_location,
               aes(x = Location, y = Rate, fill = Survival, Color = Survival)) +
  stat_summary(fun = mean, geom = "bar", position = "dodge") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 100), breaks = c(0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("green4", "red3"), 0.5)) +
  labs(title = "Survival Rate by Growing Season and Location") +
  xlab("") +
  ylab("Survival rate (%)") +
  facet_wrap( ~ Year) +
  guides(fill = "none") +
  theme_classic()

plotly::ggplotly(Sur3)
Sur4 <- ggplot(hane_survival_year_location_family,
               aes(x = Family, y = Rate, fill = Survival, color = Survival)) +
  stat_summary(fun = mean, geom = "bar", position = "dodge") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 100), breaks = c(0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("green4", "red3"), 0.5)) +
  labs(title = "Survival Rate by Growing Season, Location, and Family") +
  xlab("") +
  ylab("Survival rate (%)") +
  facet_wrap( ~ ID, ncol = 3) +
  guides(fill = "none") +
  theme_classic()

plotly::ggplotly(Sur4)

GLD

Fig1 <- ggplot(hane_gld,
                 aes(x = Family, y = GLD, fill = Family, color = Family)) +
  stat_summary(fun.data = mean_sd, geom = "errorbar", width = 0.2, color = "grey") +
  stat_summary(fun = mean, geom = "bar") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 60), breaks = c(0, 10, 20, 30, 40, 50, 60), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("#E69F00", "#56B4E9", "#009E73",
                                     "#F0E442", "#D55E00", "#CC79A7"), 0.5)) +
  labs(title = "GLD: Family Effect") +
  xlab("Family") +
  ylab("GLD (mm)") +
  annotate("text",
           x = c(1, 2, 3, 4, 5, 6),
           y = c(50),
           label = c("b", "ab", "a", "a", "ab", "ab")) +
  guides(fill = "none") +
  theme_classic()

plotly::ggplotly(Fig1)
Fig2 <- ggplot(hane_gld,
               aes(x = Year, y = GLD, fill = Location, color = Location)) +
  stat_summary(fun.data = mean_sd, geom = "errorbar", position = position_dodge(0.9), width = 0.2, color = "grey") +
  stat_summary(fun = mean, geom = "bar", position = "dodge") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 70), breaks = c(0, 10, 20, 30, 40, 50, 60, 70), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("#000000", "#999999"), 0.5)) +
  labs(title = "GLD: Year * Location") +
  xlab("Growing Season") +
  ylab("GLD (mm)") +
  annotate("text",
           x = c(1.78, 2.22, 2.78, 3.22),
           y = c(20, 20, 55, 50),
           label = c("b", "a", "a", "b")) +
  guides(fill = "none") +
  theme_classic()

plotly::ggplotly(Fig2)
Fig3 <- ggplot(hane_gld,
               aes(x = Year, y = GLD, fill = Family, color = Family)) +
  stat_summary(fun.data = mean_sd, geom = "errorbar", position = position_dodge(0.9), width = 0.2, color = "grey") +
  stat_summary(fun = mean, geom = "bar", position = "dodge") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 80), breaks = c(0, 10, 20, 30, 40, 50, 60, 70, 80), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("#E69F00", "#56B4E9", "#009E73",
                               "#F0E442", "#D55E00", "#CC79A7"), 0.5)) +
  labs(title = "GLD: Year * Family") +
  xlab("Growing Season") +
  ylab("GLD (mm)") +
  annotate("text",
           x = c(2.62, 2.78, 2.93, 3.08, 3.22, 3.38),
           y = c(45, 50, 55, 60, 55, 55),
           label = c("b", "b", "a", "a", "a", "a")) +
  guides(fill = "none") +
  theme_classic()

plotly::ggplotly(Fig3)

Height

Fig4 <- ggplot(hane_ht,
               aes(x = Family, y = Height, fill = Family, color = Family)) +
  stat_summary(fun.data = mean_sd, geom = "errorbar", width = 0.2, color = "grey") +
  stat_summary(fun = mean, geom = "bar") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 12), breaks = c(0, 2, 4, 6, 8, 10, 12), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("#E69F00", "#56B4E9", "#009E73",
                                     "#F0E442", "#D55E00", "#CC79A7"), 0.5)) +
  labs(title = "Height: Year * Location") +
  xlab("Family") +
  ylab("Height (m)") +
  annotate("text",
           x = c(1, 2, 3, 4, 5, 6),
           y = c(10),
           label = c("d", "cd", "a", "a", "bc", "ab")) +
  guides(fill = "none") +
  theme_classic()

plotly::ggplotly(Fig4)
Fig5 <- ggplot(hane_ht,
                  aes(x = Year, y = Height, fill = Location, color = Location)) +
  stat_summary(fun.data = mean_sd, geom = "errorbar", position = position_dodge(0.9), width = 0.2, color = "grey") +
  stat_summary(fun = mean, geom = "bar", position = "dodge") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 12), breaks = c(0, 2, 4, 6, 8, 10, 12), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("#000000", "#999999"), 0.5)) +
  labs(title = "Height: Year * Location" ) +
  xlab("Growing Season") +
  ylab("Height (m)") +
  annotate("text",
           x = c(3.78, 4.22, 4.78, 5.22, 5.78, 6.22, 6.78, 7.22),
           y = c(5, 4, 7, 5.5, 8.5, 7.5, 11, 7.5),
           label = c("a", "b", "a", "b", "a", "b", "a", "b")) +
  theme_classic()

plotly::ggplotly(Fig5)
Fig6 <- ggplot(hane_ht,
               aes(x = Year, y = Height, fill = Family, color = Family)) +
  stat_summary(fun.data = mean_sd, geom = "errorbar", position = position_dodge(0.9), width = 0.2, color = "grey") +
  stat_summary(fun = mean, geom = "bar", position = "dodge") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 12), breaks = c(0, 2, 4, 6, 8, 10, 12), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("#E69F00", "#56B4E9", "#009E73",
                                     "#F0E442", "#D55E00", "#CC79A7"), 0.5)) +
  labs(title = "Height: Year * Family") +
  xlab("Growing Season") +
  ylab("Height (m)") +
  annotate("text",
           x = c(2.63, 2.77, 2.92, 3.06, 3.22, 3.38,
                 3.63, 3.77, 3.92, 4.06, 4.22, 4.38,
                 4.63, 4.77, 4.92, 5.06, 5.22, 5.38,
                 5.63, 5.77, 5.92, 6.06, 6.22, 6.38,
                 6.63, 6.77, 6.92, 7.06, 7.22, 7.38,
                 7.63, 7.77, 7.92, 8.06, 8.22, 8.38
                 ),
           y = c(2.2, 2.5, 3, 3.3, 2.7, 2.9,
                 3.1, 3.3, 4.3, 4.5, 3.7, 3.9,
                 5, 5.5, 6.5, 6.8, 5.9, 6.2,
                 6.5, 7, 8.1, 8.4, 7.4, 7.8,
                 7.7, 8.5, 10.4, 10.4, 9.8, 10.2,
                 9.3, 9.2, 10.5, 10.4, 10, 10.3
                 ),
           label = c("c", "bc", "a", "a", "abc", "ab",
                     "c", "bc", "a", "a", "bc", "ab",
                     "d", "cd", "ab", "a", "bc", "ab",
                     "c", "bc", "a", "a", "b", "a",
                     "c", "c", "a", "a", "b", "ab",
                     "dc", "d", "a", "a", "bc", "ab"
                     )) +
  theme_classic()

plotly::ggplotly(Fig6)

DBH

Fig7 <- ggplot(hane_dbh,
               aes(x = Family, y = DBH, fill = Family, color = Family)) +
  stat_summary(fun.data = mean_sd, geom = "errorbar", width = 0.2, color = "grey") +
  stat_summary(fun = mean, geom = "bar") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 20), breaks = c(0, 5, 10, 15, 20), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("#E69F00", "#56B4E9", "#009E73",
                                     "#F0E442", "#D55E00", "#CC79A7"), 0.5)) +
  labs(title = "DBH: Family Effect") +
  xlab("Family") +
  ylab("DBH (cm)") +
  annotate("text",
           x = c(1, 2, 3, 4, 5, 6),
           y = c(15),
           label = c("d", "cd", "b", "b", "cb", "a")) +
  theme_classic()

plotly::ggplotly(Fig7)
Fig8 <- ggplot(hane_dbh,
               aes(x = Year, y = DBH, fill = Location, color = Location)) +
  stat_summary(fun.data = mean_sd, geom = "errorbar", position = position_dodge(0.9), width = 0.2, color = "grey") +
  stat_summary(fun = mean, geom = "bar", position = "dodge") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 20), breaks = c(0, 5, 10, 15, 20), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("#000000", "#999999"), 0.5)) +
  labs(title = "DBH: Year * Location" ) +
  xlab("Growing Season") +
  ylab("DBH (cm)") +
  annotate("text",
           x = c(0.78, 1.22, 1.78, 2.22, 2.78, 3.22, 3.78, 4.22, 4.78, 5.22),
           y = c(5, 4, 8, 7, 11, 10, 13, 12, 14, 13),
           label = c("a", "b", "a", "b", "a", "b", "a", "b", "a", "b")) +
  guides(fill = "none") +
  theme_classic()

plotly::ggplotly(Fig8)
Fig9 <- ggplot(hane_dbh,
               aes(x = Year, y = DBH, fill = Family, color = Family)) +
  stat_summary(fun.data = mean_sd, geom = "errorbar", position = position_dodge(0.9), width = 0.2, color = "grey") +
  stat_summary(fun = mean, geom = "bar", position = "dodge") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 20), breaks = c(0, 5, 10, 15, 20), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("#E69F00", "#56B4E9", "#009E73",
                                     "#F0E442", "#D55E00", "#CC79A7"), 0.5)) +
  labs(title = "DBH: Year * Family") +
  xlab("Growing Season") +
  ylab("DBH (cm)") +
  annotate("text",
           x = c(0.63, 0.77, 0.92, 1.06, 1.22, 1.38,
                 1.63, 1.77, 1.92, 2.06, 2.22, 2.38,
                 2.63, 2.77, 2.92, 3.06, 3.22, 3.38,
                 3.63, 3.77, 3.92, 4.06, 4.22, 4.38,
                 4.63, 4.77, 4.92, 5.06, 5.22, 5.38
           ),
           y = c(3.2, 3.5, 5, 5.3, 4.7, 4.9,
                 6.1, 6.3, 7.8, 7.9, 7.5, 8.5,
                 8.9, 9.3, 10.1, 10.2, 10.8, 11.8,
                 11.1, 11.3, 11.9, 12.3, 12.8, 13.2,
                 12.7, 12.7, 12.9, 13.1, 13.5, 14.7
           ),
           label = c("c", "c", "ab", "a", "bc", "a",
                     "d", "cd", "b", "b", "bc", "a",
                     "c", "c", "b", "b", "b", "a",
                     "d", "cd", "bcd", "bc", "b", "a",
                     "b", "b", "b", "b", "b", "a"
           )) +
  guides(fill = "none") +
  theme_classic()

plotly::ggplotly(Fig9)

Crown

Fig10 <- ggplot(hane_crown,
                aes(x = Family, y = Crown, fill = Family, color = Family)) +
  stat_summary(fun.data = mean_sd, geom = "errorbar", width = 0.2, color = "grey") +
  stat_summary(fun = mean, geom = "bar") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 10), breaks = c(0, 2, 4, 6, 8, 10), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("#E69F00", "#56B4E9", "#009E73",
                                     "#F0E442", "#D55E00", "#CC79A7"), 0.5)) +
  labs(title = "Crown Height: Family Effect") +
  xlab("Family") +
  ylab("Crown Height (m)") +
  annotate("text",
           x = c(1, 2, 3, 4, 5, 6),
           y = c(8),
           label = c("b", "b", "a", "a", "a", "a")) +
  guides(fill = "none") +
  theme_classic()

plotly::ggplotly(Fig10)
Fig11 <- ggplot(hane_crown,
               aes(x = Year, y = Crown, fill = Location, color = Location)) +
  stat_summary(fun.data = mean_sd, geom = "errorbar", position = position_dodge(0.9), width = 0.2, color = "grey") +
  stat_summary(fun = mean, geom = "bar", position = "dodge") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 10), breaks = c(0, 2, 4, 6, 8, 10), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("#000000", "#999999"), 0.5)) +
  labs(title = "Crown Height: Year * Location" ) +
  xlab("Growing Season") +
  ylab("Crown Height (m)") +
  annotate("text",
           x = c(0.78, 1.22, 1.78, 2.22, 2.78, 3.22),
           y = c(3.2, 1.7, 6.5, 4, 7.3, 5),
           label = c("a", "b", "a", "b", "a", "b")) +
  guides(fill = "none") +
  theme_classic()

plotly::ggplotly(Fig11)
Fig12 <- ggplot(hane_crown,
               aes(x = Year, y = Crown, fill = Family, color = Family)) +
  stat_summary(fun.data = mean_sd, geom = "errorbar", position = position_dodge(0.9), width = 0.2, color = "grey") +
  stat_summary(fun = mean, geom = "bar", position = "dodge") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 10), breaks = c(0, 2, 4, 6, 8, 10), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("#E69F00", "#56B4E9", "#009E73",
                                     "#F0E442", "#D55E00", "#CC79A7"), 0.5)) +
  labs(title = "Crown Height: Year * Family") +
  xlab("Growing Season") +
  ylab("Crown Height (m)") +
  annotate("text",
           x = c(1.63, 1.77, 1.92, 2.06, 2.22, 2.36,
                 2.63, 2.77, 2.92, 3.06, 3.22, 3.36
           ),
           y = c(2.5, 3.4, 6, 6.5, 5.5, 5.7,
                 4, 3.9, 7.2, 7.3, 6.1, 6.6
           ),
           label = c("b", "b", "a", "a", "a", "a",
                     "b", "b", "a", "a", "a", "a"
           )) +
  guides(fill = "none") +
  theme_classic()

plotly::ggplotly(Fig12)